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ABSTRACT 


Modern  naval  warfare  has  been  increasingly  dependent  upon  the 
acoustic  silencing  of  the  participants.  Constrained  viscoelastic  layer 
damping  of  vibrating  elements  is  one  method  which  can  be  used  to  meet 
acoustic  silencing  goals.  This  paper  considers  constrained  viscoelastic 
layer  damping  treatments  applied  to  a  thick  aluminum  plate,  including 
single  layer,  double  layer,  a  milled  pocket  plate,  and  a  milled  “floating 
element”  configuration.  High  modal  damping  values  were  obtained  for 
each  damping  configuration.  The  Modal  Strain  Energy  method,  using 
finite  element  analysis  to  estimate  modal  loss  factors,  was  investigated  for 
use  as  a  tool  in  constrained  viscoelastic  layer  damping  design.  A 
comparison  of  experimentally  measured  frequency  response  and  modal 
loss  factors  with  those  predicted  by  the  modal  strain  energy  method  is 
presented  to  confirm  the  possible  use  of  the  modal  strain  energy  method  as 
a  design  tool. 
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L  INTRODUCTION 


A.  BACKGROUND 

Modem  naval  warfare,  especially  undersea  warfare,  depends  heavily 
on  the  vessel  being  acoustically  silent.  A  major  source  of  radiated  noise  is 
the  vibration  of  shipboard  components.  The  reduction  of  these  vibrations  is 
of  utmost  importance  if  a  ship  is  to  accomplish  its  mission.  One  method  of 
vibration  damping  that  shows  promise  in  damping  over  a  broad  spectrum 
of  low  frequency  vibration  is  constrained  viscoelastic  layer  damping.  The 
constrained  viscoelastic  layer  method  uses  shear  deformation  in 
viscoelastic  material  to  absorb  and  dissipate  the  vibrational  energy  of  the 
system  in  question.  Unfortunately  the  design  and  analysis  of  such 
constrained  viscoelastic  layer  systems  is  difficult,  due  in  part  to  the 
following: 

•  The  material  properties  of  viscoelastic  damping  materials  vary 
greatly  with  temperature  and  frequency. 

•  Closed  form  solutions  to  the  equations  of  motion  exist  only  for  beams 
and  plates  with  simple  boundary  conditions. 

•  The  necessity  for  large  amounts  of  computer  storage  and  CPU  time  to 
conduct  finite  element  analysis  of  constrained  viscoelastic  layer 
damping  systems  because  of  a  large  number  of  elements  and  the 
variation  of  material  properties  with  frequency. 
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A  finite  element  technique  developed  by  Johnson  and  Kienholz 
known  as  the  Modal  Strain  Energy  (MSE)  method  uses  structural  strain 
energy  to  approximate  the  damping  of  a  structure  with  an  applied 
constrained  viscoelastic  damping  system  [Ref.  1].  Previous  work  by  Maurer 
addressed  the  effectiveness  of  the  MSE  method  for  two  damped  plate 
configurations:  1)  a  simple  sandwich  configuration,  and  2)  a  plate  with  a 
milled  pocket  with  damping  material  inserted  and  a  welded  cover  plate 
acting  as  a  constraining  layer  [Ref. 2].  Difficulties  with  the  welded  cover 
plate  were  reported  because  during  welding,  the  cover  plate  warped  and 
delaminated  itself  from  the  damping  material,  resulting  in  negligible 
damping  [Ref.  2], 

B.  PURPOSE 

For  certain  naval  applications  the  components  to  be  damped  will  be 
thick  in  construction  and  may  be  exposed  to  an  unfriendly  environment. 
This  paper  addresses  the  experimental  testing  and  analysis  of  four  thick 
aluminum  plates,  each  with  a  different  constrained  viscoelastic  layer 
damping  treatment. 

This  research  further  investigated  the  pocket  plate  configuration  and 
also  investigated  a  constrained  viscoelastic  layer  damping  system  using 
two  separate  layers  of  damping  material.  In  addition,  a  second  milled  plate 
using  a  “floating  element”  in  conjunction  with  constrained  layer  damping 
was  tested.  The  MSE  method  is  used  to  analyze  each  of  the  structures,  and 
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its  accuracy  and  usefulness  as  a  possible  design  tool  are  investigated. 
Although  the  intended  use  of  the  MSE  method  is  the  design  of  complex 
damping  structures,  simple  plate  geometries  were  used  to  facilitate  the 
experimental  and  computational  effort. 
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II.  THEORY 


A.  VISCOELASTIC  MATERIAL 

Viscoelastic  materials  of  interest  for  general  naval  applications  are 
polymeric  compounds  made  up  of  long  molecular  chains.  These  molecular 
chains  can  be  strongly,  or  weakly,  linked  together,  depending  on  their 
chemical  composition  and  processing.  The  damping  characteristics  of 
viscoelastic  arise  from  the  deformation  and  recovery  of  the  polymer 
network.  Material  properties  of  a  viscoelastic  material  vary  with 
temperature  and  frequency.  As  such,  the  damping  characteristics  of  a 
system  will  vary  as  its  operating  environment  changes.  [Ref.  3] 

Temperature  will  have  the  greatest  effect  on  the  material  properties 
of  damping  materials  [Ref.  3].  This  effect  is  shown  in  Figure  2.1,  where 
four  distinct  regions  are  observed.  The  lowest  temperature  region  is  the 
glassy  region  where  the  material’s  storage  modulus  is  at  its  maximum 
value,  and  the  loss  factor  is  at  a  minimum.  In  the  glassy  region  the 
modulus  decreases  slowly  with  temperature  increase,  whereas  the  loss 
factor  increases  rapidly  with  temperature.  The  second  region  is  the 
transition  region  where  the  modulus  decreases  rapidly  with  increasing 
temperature  and  the  loss  factor  reaches  its  peak  value.  The  third  region  is 
the  rubbery  region  where  both  the  modulus  and  loss  factor  are  at  low  values 
and  show  little  variation  with  temperature.  The  fourth,  and  last,  region  is 
the  flow  region  and  characterizes  the  behavior  of  some  materials,  mostly 
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ceramics,  at  high  temperatures.  It  should  be  noted  that  the  transition 
region  may  vary  in  width  from  20  °C  up  to  a  width  of  200  °C.  [Ref.  3] 

The  effect  of  frequency  on  viscoelastic  materials  is  not  as  great  as  that 
of  temperature.  The  modulus  of  the  viscoelastic  always  increases  with 
increasing  frequency.  The  loss  factor  will  initially  increase  with  frequency, 
then  peak,  and  subsequently  decrease  as  frequency  increases.  A  plot  of 
storage  modulus  and  loss  factor  versus  frequency  is  shown  in  Figure  2.2.  It 
should  be  noted  that  this  plot  is  over  a  range  of  approximately  ten  decades, 
and  hence  it  becomes  obvious  that  a  temperature  change  of  a  couple  degrees 
will  have  a  much  greater  effect  on  damping  than  a  minor  change  in 
frequency.  [Ref.  3] 


Figure  2.1.  Variation  of  viscoelastic  material 
properties  with  temperature  [Ref.  3]. 


5 


Figure  2.2.  Variation  of  viscoelastic  material  properties 
with  frequency  and  constant  temperature  [Ref.  3]. 


Linear  viscoelastic  materials  behave  in  a  hysteretic  manner  under 
cyclic  excitation.  Therefore,  the  mechanical  behavior  of  a  viscoelastic 
material  during  steady  state  vibration  is  best  described  by  using  a  complex 
stiffness,  k*  [Ref.  3]. 

k*  =  k(l  +  irj)  (2.1) 

where, 

r|  =  material  loss  factor 


The  use  of  a  complex  stiffness  then  leads  to  the  use  of  a  complex  Young’s 
modulus  and  shear  modulus  [Ref.  3]. 

E*  =  E(1  +  irt)  (2.2) 

G*  =  Gd  +  ir|)  (2.3) 

This  concept  of  the  complex  modulus  is  used  in  subsequent  analysis. 

Viscoelastic  material  properties  are  commonly  displayed  using  a 
“reduced  frequency  nomogram.”  The  reduced  frequency  nomogram 
displays  the  variation  of  the  viscoelastic  material’s  loss  factor  and  modulus 
with  temperature  and  frequency.  The  “reduced  frequency”,  f«t>  is  an 

empirically  determined  function  that  accounts  for  the  viscoelastic’s 
temperature  and  frequency  dependence,  and  allows  data  for  wide  range  of 
temperature  and  frequencies  to  be  plotted  on  the  same  graph  [Ref  4].  The 
reduced  frequency  nomogram  for  3M  ISD  -  112  is  shown  in  Figure  2.3.  To 
find  the  loss  factor  and  modulus  using  the  nomogram,  enter  with  the 
desired  temperature  and  frequency.  Follow  the  frequency  line  horizontally 
and  the  temperature  line  diagonally  down  the  page  until  the  two  intersect. 
Then  go  vertically  up  or  down  to  intersect  the  shear  modulus  or  loss  factor 
curves.  Finally,  read  the  value  of  the  shear  modulus  or  loss  factor 
horizontally  from  the  scale  on  the  left  [Ref.  4]. 
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DAMPING  PROPERTIES 
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R  CONSTRAINED  VISCOELASTIC  LAYER  DAMPING 

A  simple  constrained  layer  damping  treatment  consists  of  a  base 
layer  (the  structure  to  be  damped),  a  damping  layer,  and  the  constraining 
layer.  This  configuration  is  shown  in  Figure  2.4  with  the  thicknesses  of  the 
damping  and  constraining  layers  exaggerated  for  clarity. 


Figure  2.4.  Simple  constrained  viscoelastic  layer  configuration. 

The  physical  mechanism  of  damping  can  be  explained  by  referring  to 
Figure  2.4.  When  the  base  layer  is  deformed  in  a  mode  of  vibration,  the 
surface  away  from  the  neutral  axis  elongates,  stretching  the  viscoelastic 
material.  The  top  layer,  being  a  stiff  elastic  material,  tends  not  to  elongate, 
and  thereby  “constrains”  the  viscoelastic  material.  Consequently,  the  cyclic 
motions  of  vibration  induce  a  cyclic  shearing  strain  in  the  viscoelastic. 
This  cyclic  shearing  strain,  together  with  its  associated  hysteresis  loop 
cause  the  vibrational  energy  to  be  dissipated  as  heat.  For  the  constraining 
layer  to  be  effective,  its  stiffness  should  not  exceed  that  of  the  base  layer. 
[Ref.  4  &  5] 
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C.  SYSTEM  EQUATIONS  OF  MOTION 


Continuous  systems,  such  as  plates,  possess  distributed 
characteristics  of  mass,  damping,  and  stiffness.  Classical  vibration 
analysis  of  such  systems  involves  the  formation  of  a  mathematical  model 
that  discretizes  the  system  into  a  finite  number  of  components  in  order  to 
approximate  the  total  system.  Such  a  formulation  results  in  the  following 
equation: 

[M]  (x(t)}  +  [C]  (x(t)}  +  [K]  { x(t) }  =  { F(t) }  (2.4) 

where, 

[M]  =  system  mass  matrix 
[C]  =  system  damping  matrix 
[K]  =  system  stiffness  matrix 
(F(t)}  =  external  excitation  vector 
(x(t)}  =  displacement  vector 

For  an  undamped  system  without  excitation,  the  above  equation  reduces  to 
the  eigenvalue  problem. 


[M]{x(t)}+[K]{x(t)}  =  0  (2.5) 

This  equation  is  then  transformed  to  modal  space  using  the  linear 
transformation: 


(x(t)}  =[d>]  (q(t)j 


(2.6) 


where, 


[<{>]  =  modal  matrix 

(q(t)}  =  modal  response  vector 

Using  this  linear  transformation  the  equation  of  motion  can  then  be  solved 
for  the  undamped  modal  frequencies  and  mode  shapes. 

To  solve  for  the  frequency  response  of  a  damped  system,  the  linear 
transformation  is  applied  to  equation  (2.4): 

[M]  [<t>]  |q(t))  +  [C]  [<>]  {q(t»  +  [K]  [<J>]  {q(t)>  =  (F(t)}  (2.7) 

Assuming  that  the  damping  matrix  [C]  is  proportional  to  a  linear 
combination  of  the  stiffness  matrix  [K]  and  mass  matrix  [M],  the  damping 
matrix  can  then  be  diagonalized  using  the  same  linear  transformation 
used  to  diagonalize  [K]  and  [M]  in  equation  (2.5)  above.  The  diagonal  terms 
of  the  damping  matrix  then  become  (rjj.coj),  where  Tii  equals  the  modal  loss 
factor  and  is  the  natural  frequency  of  the  i^1  mode  [Ref.  1]  .  Using  this 
approximate  diagonal  damping  matrix  results  in  a  system  of  uncoupled 
modal  equations  of  motion: 

cji(t)  +  riicq<j(t)  +  (opqi(t)  =  fi(t)  (2.8) 


where, 

ijj(t)  =  modal  acceleration  of  i^  mode 
q(t)  =  modal  velocity  of  i^1  mode 
qj(t)  =  modal  displacement  of  i^1  mode 
"Hi  =  modal  loss  factor  of  i^  mode 


coj  =  natural  frequency 
fi(t)  =  modal  force  in  mode 
j  =  ^ 

Assuming  that  a  sinusoidal  excitation  produces  a  sinusoidal  response, 

{fU)}  =  (f!ei“*  {q(t)}  =  {Q}e><*  (2.9) 

the  response  for  the  i^  mode  is  then  solved  to  be: 

q.=  — -  (2.10) 

cof  -  or  +  jcDiTtito 

Subsequently,  the  response  of  the  physical  system  can  be  found  using: 

(x(t)}  =  [<J>]  {Q)e>°*  (2.11) 

D.  MODAL  STRAIN  ENERGY  METHOD 

The  equations  of  motion  used  to  define  the  response  of  a  system  with 
viscoelastic  materials  need  a  complex  eigenvalue  analysis.  However,  the 
actual  solution  of  these  equations  may  be  quite  difficult.  This  is  especially 
true  when  the  system  to  be  analyzed  is  comprised  of  materials  whose 
properties  vary  with  both  temperature  and  frequency.  Finite  element 
techniques  are  generally  used  to  compute  the  response  of  complicated 
systems.  However,  for  the  case  of  varying  material  properties  many  time 
consuming  and  costly  runs  must  be  made  with  the  material  properties 


changing  at  each  frequency  increment  [Ref.  1].  In  addition  to  the  costly 
analysis  of  a  single  design  configuration,  changes  in  design  options,  design 
requirements,  of  the  search  for  an  optimum  design  can  make  the  expense 
of  finite  element  analysis  too  great.  The  development  of  the  Modal  Strain 
Energy  (MSE)  method  by  Johnson  and  Kienholz,  however,  makes  the  finite 
element  analysis  of  complex  viscoelastically  damped  structures  a  viable 
option  [Ref.  1]. 

The  MSE  method  assumes  that  a  damped  structure  can  be 
represented  by  the  normal  modes  of  the  associated  undamped  system  if 
appropriate  damping  terms  are  inserted  into  the  uncoupled  modal 
equations  of  motion  (eqn  2.8)  [Ref.  1].  The  MSE  method  further  assumes 
that  the  damping  matrix  of  equation  (2.4)  can  be  diagonalized  by  the  same 
real  modal  matrix  that  diagonalizes  the  system  mass  and  stiffness 
matrices  of  equation  (2.4). 

The  modal  loss  factors  of  equation  (2.8)  are  calculated  using  the 
undamped  mode  shapes  and  material  loss  factors  for  each  material.  Since 
the  material  loss  factors  of  the  structure  to  be  damped  and  the  constraining 
layers  are  quite  small  compared  to  those  of  common  viscoelastic  cores,  the 
modal  loss  factors  can  be  estimated  using  the  following  equation  [Ref.  1]: 

Tl(r)  =  Tl(^  (2.12) 


where, 

Tl(r)  =  system  loss  factor  in  the  r^  mode 


Ti(vr)  =  material  loss  factor  of  viscoelastic  core  at  the  r^ 
natural  frequency 

Vyr)  =  strain  energy  in  the  viscoelastic  core  at  r^  mode 
V(r)  =  strain  energy  in  the  entire  structure  at  r^  mode 

The  modal  strain  energies  can  be  obtained  for  finite  element  analysis,  and 
are  a  standard  output  option  of  the  NASTRAN  finite  element  code  [Ref.  7], 
The  modal  frequency  response  of  the  structure  is  then  calculated 
using  the  modal  loss  factors  found  in  equation  (2.12).  When  computing  the 
modal  frequency  response  of  a  damped  structure,  the  modal  properties  in 
the  system  matrices  are  assumed  to  be  constant.  However,  viscoelastic 
materials  have  storage  moduli  which  are  frequency  dependent.  To  account 
for  this  frequency  dependence,  Johnson  and  Kienholz  devised  the  following 
correction  factor  to  be  applied  to  the  modal  loss  factors  calculated  in 
equation  (2.12)  [Ref.  8]. 


(2.13) 


where, 

Tt(cr)  =  corrected  modal  loss  factor  at  the  r^1  mode 
G2(fr)  -  viscoelastic  shear  modulus  at  the  r^  modal 
frequency 

G2,ref  =  reference  viscoelastic  shear  modulus  used  in  the 
frequency  response  calculation 


IIL  DESIGN  OF  DAMPED  PLATES 


A .  GENERAL  SPECIMEN  CONFIGURATIONS 

For  experimental  testing  and  analysis  purposes,  four  different,  yet 
related,  constrained  layer  damping  treatments  were  selected  in  addition  to 
an  undamped  “reference”  plate.  Two  of  the  damping  configurations  were 
simple  sandwich  treatments  consisting  of  one  and  two  viscoelastic  layers 
respectively.  Another  treatment  was  the  “pocket  plate”  which  was 
previously  investigated  by  Maurer  [Ref.  2].  The  pocket  plate  was  made  from 
a  solid  plate  which  was  then  milled  to  accept  damping  material  and  a  cover 
plate.  The  final  damping  treatment  was  a  “floating  element”  configuration 
consisting  of  a  solid  plate  milled  to  accept  a  double  layer  of  damping 
material  and  a  welded  cover  plate.  Section  views  of  these  damping 
configurations  are  shown  in  Figure  3.1. 

The  purpose  of  the  pocket  plate  is  to  protect  the  viscoelastic  material 
from  materials  such  as  oil  and  salt  water,  which  may  harm  the 
viscoelastic.  Since  previous  attempts  at  using  a  welded  cover  plate  were 
unsuccessful  due  to  the  heat  of  welding  causing  a  delamination  between 
the  damping  material  and  cover  plate  [Ref.  2],  it  was  decided  to  move  the 
viscoelastic  material  from  the  welding  point  by  recessing  it  into  a  shallow 
pocket  of  its  own  as  shown  in  Figure  3.1. 

The  “floating  element”  concept  evolved  because  the  welded  cover  plate 
of  the  pocket  plate  configuration  does  not  produce  the  damping  reaction  that 
a  true  constraining  layer  would  provide.  If  the  cover  plate  is  welded  to  the 


surrounding  structure  it  cannot  deform  in  bending  as  much  as  an 
unwelded  constraining  layer,  thereby  causing  a  reduction  in  the  damping 
capability  of  the  system.  By  using  a  piece  of  metal  in  the  milled  pocket  with 
dimensions  slightly  smaller  than  the  surrounding  pocket,  along  with  two 
layers  of  viscoelastic  and  a  welded  cover  plate,  a  true  constraining  layer 
effect  should  be  obtained.  This  configuration  is  shown  in  Figure  3.1 . 


Single  layer  configuration 


Double  layer  configuration 


ir/r>r-?-r>/rr7-7-77c 


Milled  “pocket  plate” 


Plate  with  “floating  element” 


Figure  3.1.  Generalized  damping  treatment  configurations. 


In  order  to  approximate  a  possible  system  to  be  damped,  a  plate  with 
large  dimensions  was  selected.  The  dimensions  of  the  plates  used  for  the 
damping  treatments  are  114.3  cm  (45  in)  in  length  and  38.1  cm  (15  in)  in 
width.  In  addition,  possible  naval  applications  for  this  type  of  damping 
would  probably  consist  of  thick  plate  members.  For  this  reason  it  was  also 
decided  that  thick  plates  would  be  used  for  the  damping  treatments.  All  the 


base  plates  and  constraining  layers  were  made  of  a  standard  6061-T6 
aluminum  alloy.  The  design  and  selection  of  viscoelastic  layer,  base  layer, 
and  constraining  layer  thicknesses  is  discussed  in  the  following  sections. 

B.  DESIGN  OF  THE  SINGLE  DAMPING  LAYER  CONFIGURATION 

In  an  attempt  to  approximate  system  loss  factors  and  hence 
determine  viscoelastic  and  constraining  layer  thicknesses  for  maximum 
damping,  a  method  de /eloped  by  Nashif  [Ref.  9]  based  on  an  analysis  of 
simple  sandwich  damping  systems  by  Ross,  Kerwin,  and  Ungar  [Ref.  10] 
was  used.  The  Ross  -  Kerwin  -  Ungar  (RKU)  equations  are  base  of  the 
analysis  of  the  simple  sandwich  system  shown  in  Figure  3.2. 


CONSTRAINING  LAYER 


DAMPING  LAYER 


BASE  LAYER 


Figure  3.2.  Elements  of  a  simple  sandwich  damping  system. 


To  find  the  loss  factors  of  the  damped  system,  the  flexural  rigidity  of 
the  system  must  first  be  determined.  For  the  above  system,  the  flexural 
rigidity,  El,  is  written  as  [Ref.  9]: 


El  =  Eill  +  +  EsMl  +  E  HlD2 

12  12  12 

+  E2H2(H21  -  D)2+  E3Ha(H3i  -  D)2 


Sfp  +  ®^H21  Dl  +  EaHalHa,  D)]Hai_J> 

L  12  2  J  1  +  g 

(3.1) 

where, 

E2HJH21  -  ^L)  +  g(E2H2H2i  +  E3H3H31) 

D  =  '  2  ' 

EiHj  +  pp  +  g(E1H1  +  E2H2  +  E3H3) 

(3.2) 

Hai  =  Sl±Jil  +  h2 

(3.3) 

h21  =  h.+2H2 

(3.4) 

g=  °2  2 

e3h3h2k2 

(3.5) 

E  =  Young’s  modulus 
G  =  Shear  modulus 
I  =  second  moment  of  area 
H  =  thickness  of  member 
K2  =  wave  number 

Subscripts  refer  to  the  layers  labeled  in  Figure  3.1 
No  subscript  refers  to  the  composite  system 


For  a  simply  supported  plate  the  wave  numbers  and  modal  frequencies  are 
found  using  [Ref.9]: 


—  Knm 


EH3gc 
12(l  -  v2)Hp 


(3.6) 


where, 

a  =  semi— wave  length  of  the  plate 
b  =  semi-wave  width  of  the  plate 
v  =  Poisson’s  ratio  of  the  composite  plate 
p  =  density  of  the  composite  plate 
gc  =  gravitational  constant 


(3.7) 


To  introduce  damping  into  the  equations  it  is  necessary  to  use  the 
complex  modulus  concept  expressed  in  Section  II.  Substituting  the 
appropriate  expressions  for  the  complex  shear  and  Young’s  modulus  into 
equations  (3.1),  and  assuming  that  damping  in  the  base  layer,  (ri^),  is 

small,  and  that  the  extensional  stiffness  of  the  damping  layer  is  small 
(since  E2  «  E^  and  E2  «  E3),  the  following  expressions  can  be  arrived  at 

[Ref.9]: 


EH3  =  ExH?  +  E3Hj]  +  —12—  (a  -  p  -  3)re 

c2  +  d2 


1  9 


(3.8) 


(3.9) 


EH3!!  =  E3H3TI3  +  — 12—  (a  -  p  -  3)m 
c2  + 


where, 

a  =  gEiHxEaHsHiilcfl  -  TfeTts)  +  dft2  +  *n3)  +  j[c(rj2  +  n3)]}  (3.10) 

p  =  EiHiE^^si  [c  +  drj2  +  j(cT]2  -  d)]  (3.1 1 ) 


d  =  2gE2H2E3H3H2iH3i 


c(l  -2Tl2ri3-T}D  +  d(2Tl2  +  Tl3-Tl|Tl3)  \ 
j[c(2rj2  +  T13  -  Tiins)  -  d(l  -  2ri2ri3  -  til)]) 


(3.12) 


c  =  EiHi(l  +  g)  +  gE3H3(l  -r\2r\3)  (3.1 3) 

d  =  gE1HiTj2  +  gE3Ha(Ti2  +  ti3)  (3.14) 

j  =  (3.15) 


These  equations  were  then  applied  to  estimate  the  loss  factors  of  the  simple 
three— layer  sandwich  plate.  The  equations  can  be  simplified  by  assuming 
that  damping  in  the  constraining  layer  (TI3)  is  negligible  [Ref.9]. 

Since  the  boundary  conditions  for  the  plate  used  in  this  research 
(free-free-free-free)  do  not  correspond  to  the  simply  supported  conditions  on 
which  equation  (3.7)  is  based,  modal  frequencies  for  the  free-free  case  were 
estimated  using  results  from  finite  element  analysis  of  a  free-free  plate  and 
equation  (3.6).  The  natural  frequencies  of  an  undamped  plate  with  the 
dimensions  previously  given,  and  a  thickness  of  1.91  cm  (0.75  in)  were 
found  using  a  normal  mode  extraction  in  NASTRAN.  By  substituting  these 
modal  frequencies  into  equation  (3.6)  and  estimate  of  the  wave  parameter 
for  each  mode,  K^m,  was  obtained.  Then,  by  using  an  iterative  procedure 
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outlined  in  Reference  [9],  modal  loss  factors  were  estimated  for  different 
layer  thicknesses  over  a  temperature  range  of  0.0  °C  to  37.8  °C  (30  °F  to 
100  °F). 

The  previous  equations  are  easily  programmed  to  compute  loss 
factors  for  a  wide  variety  of  conditions.  The  variation  of  viscoelastic 
material  properties  with  temperature  and  frequency  was  accounted  for 
using  a  curve-fit  to  the  reduced  frequency  nomogram  developed  by  Drake 
[Ref.  11].  The  material  data  for  the  following  curve-fit  equations  is  from  the 
University  of  Dayton  Research  Institute  [Ref.  12]. 


2 


logi0{M)  =  logio(ML)  + 


(3.16) 


log10(ETA)  =  logio(ETAFROL) 

+  lc[A(SL  +  SH)  +  (SL -  SH/l  -  Vl  +  A2)]  ('  } 

jL 

logiotFR)  ■  logio(F)  -  (625  ^  j  <318> 


A  logi  o(FR)  -  logt  o(FROL) 

A"  C 

where, 

M  =  viscoelastic  modulus 

ETA  =  viscoelastic  material  loss  factor 

FR  =  reduced  frequency  (Hz) 

F  =  frequency  (Hz) 

T  =  temperature  (°F) 


(3.19) 


and, 


T0  =  40  °C  (104  °F) 

FROM  =  2.0x1 04  Hz 

MROM  =  4.75x1 06  Pa  (688.94  psi) 

n  =  0.275 

ML  =  6.0x1 04  Pa  (8.7  psi) 

ETAFROL  =  1.08 
SL  =  0.45 
SH  =  -0.55 
FROL  =  5000  Hz 
C  =  2.5 

In  addition  to  the  above  constants,  a  Poisson’s  ratio  of  0.49  and  a  density  of 
0.909  gram  per  cubic  centimeter  (0.035  Ibm/in®)  was  used  for  ISD-112 
[Ref.12],  The  following  material  properties  were  used  for  6061-T6 
aluminum  [Ref.13]. 

E  =  70  GPa  (10x10®  psi) 
v  =  0.33 

p  =  2.7  gm/cm®  (0.0968  lbm/in®) 

Using  the  previous  equations  and  material  properties,  a  computer 
program  was  written  to  compute  estimated  modal  loss  factors  and 
frequencies  for  a  variety  of  base  layer,  viscoelastic  layer,  and  constraining 
layer  thicknesses.  A  listing  of  this  program  appears  in  Appendix  A. 
Modal  loss  factors  were  computed  for  base  layer  thicknesses  of  9.53  mm 
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(0.375  in)  to  19.05  mm  (0.75  in)  in  3.18  mm  (0.125  in)  increments.  For  each 
base  layer  thickness,  the  viscoelastic  thickness  was  varied  from  0.38  mm 
(0.015  in)  to  1.52  mm  (0.060  in)  in  0.38  mm  increments,  and  the 
constraining  layer  thickness  was  varied  from  1.59  mm  (0.0625  in)  to 
6.35  mm  (0.25  in)  in  1.59  mm  increments.  In  addition,  loss  factors  were 
also  computed  for  a  viscoelastic  thickness  of  0.127  mm  (0.005  in).  From  the 
results  of  the  analysis,  a  carpet  plot  [Ref.  3],  was  made  for  each  of  the  base 
layer  conditions.  The  carpet  plot  reflects,  for  the  the  first  mode,  the 
maximum  loss  factor  and  its  corresponding  temperature  for  each 
viscoelastic  layer/constraining  layer  thickness  configuration.  The  carpet 
plot  for  a  base  layer  thickness  of  12.7  mm  (0.50  in)  is  shown  in  Figure  3.3. 
Based  on  the  carpet  plots  and  a  desire  for  maximum  damping,  as  well  as  a 
system  which  could  be  moved  easily,  a  base  layer  thickness  of  12.7  mm 
(0.50  in),  a  viscoelastic  thickness  of  0.38  mm  (0.015  in),  and  a  constraining 
layer  thickness  of  6.35  mm  (0.25  in)  was  selected.  The  total  system 
thickness  was  approximately  19.05  mm  (0.75  in).  This  total  system 
thickness  would  be  maintained  for  all  subsequent  damping  configurations. 

To  maintain  continuity  between  damping  systems,  the  milled  “pocket 
plate”  was  given  a  viscoelastic  thickness  of  0.38  mm  and  a  cover  plate 
thickness  of  6.35  mm  for  a  total  system  thickness  of  19.05  mm  (0.75  in).  In 
order  to  keep  the  heat  of  welding  away  from  the  viscoelastic  material,  the 
ISD-112  was  recessed  into  a  shallow  pocket  as  indicated  previously  in 
Figure  3.1 ,  and  as  shown  in  the  pocket  plate  system  arrangement  in 
Figure  3.4.  Detail  drawings  of  the  pocket  plate  are  shown  in  Appendix  C. 
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Maximum  System  Loss  Factor 


Figure  3.3.  Carpet  plot  for  a  base  layer  thickness,  HI  =  12.7  mm. 


C.  DESIGN  OF  THE  DOUBLE  DAMPING  LAYER 


The  design  of  the  double  layer  damping  system  was  accomplished 
using  a  modification  of  the  RKU  analysis  used  in  the  previous  section.  The 
RKU  equations  are  used  by  working  from  the  top  layer  of  the  damping 
system  down  towards  the  base  layer.  As  shown  in  Figure  3.5,  the  H3' 
constraining  layer  along  with  the  H2'  viscoelastic  layer  are  combined  with 
the  HI'  layer  to  form  a  three-layer  system.  Using  the  RKU  equations,  the 
stiffness  of  this  system  is  computed  and  considered  to  be  the  equivalent 
stiffness  of  the  top  three  layers  of  the  total  constrained  layer  damping 
system.  The  top  three  layers  were  then  considered  as  a  single  layer  with 
the  equivalent  stiffness  previously  calculated,  and  the  RKU  equations  were 
again  applied  to  compute  estimated  modal  loss  factors  for  the  entire  double 
layer  damping  system.  [Ref.3] 


CONSTRAINING  LAYER 
VISCOELASTIC 


CONSTRAINING  LAYER 
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H2’ 
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HI 


Figure  3.5.  General  configuration  of  the  double 
constrained  layer  system. 
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To  maintain  continuity  among  all  the  damping  configurations,  a 
base  layer  thickness  of  12.77  mm  (0.50  in)  was  chosen,  and  a  total  system 
thickness  of  19.05  mm  (0.75  in)  was  maintained.  A  design  for  high 
damping  was  then  selected  by  computing  modal  loss  factors  for  the 
constraining  and  damping  layer  thickness  combinations  shown  in  Table 
3.1 .  The  estimated  modal  loss  factors  for  the  first  mode  of  vibration  in  each 
configuration  are  plotted  as  shown  in  Figure  3.6.  A  listing  of  the  program 
used  to  compute  the  loss  factors  is  in  Appendix  B. 

From  the  data  presented  in  Figure  3.6,  viscoelastic  thickness  of 
0.38  mm  (0.015  in)  and  constraining  layer  thicknesses  of  3.18  mm  (0.125  in) 
were  selected  for  the  double  layer  configuration.  This  particular 
configuration  estimates  high  damping  over  a  wider  temperature  range 
than  the  other  thickness  combinations. 


TABLE  3  X  THICKNESSES  USED  IN  CALCULATION  OF  DOUBLE 
LAYER  MODAL  LOSS  FACTORS. 


Confi  miration 

1 

2 

2 

4 

5 

Q 

H2  (mm) 

0.38 

0.38 

0.38 

0.76 

0.38 

1.14 

HI'  (mm) 

3.18 

2.38 

3.18 

2.38 

2.38 

2.38 

H2'  (mm) 

0.38 

0.38 

0.38 

0.76 

1.14 

0.38 

H3'  (mm) 

3.18 

3.18 

2.38 

2.38 

2.38 

2.38 
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Figure  3.6.  Modal  loss  factors  for  double  constrained  layer 
damping  configuration. 

The  milled  “floating  element”  plate  uses  the  same  viscoelastic  and 
constraining  layer  thicknesses  as  the  simple  two-layer  configuration.  In  a 
design  similar  to  that  of  the  pocket  plate,  the  floating  element  and  both 
layers  of  viscoelastic  are  recessed  into  a  milled  opening  as  shown 
previously  in  Figure  3.1  and  further  described  in  the  floating  element 
system  configuration  of  Figure  3.7. 
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Figure  3.7.  Floating  element  system  configuration. 
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IV.  EXPERIMENTAL  RESULTS 


A.  TESTING  ARRANGEMENT 

Experimental  testing  was  performed  on  each  of  the  four  damping 
configurations  and  the  undamped  reference  plate.  In  order  to  approximate 
the  free-free-free-free  boundary  condition,  each  plate  was  suspended  from 
the  roof  of  the  testing  chamber  using  elastic  cords  as  shown  in  Figure  4.1 . 
All  of  the  tests  were  performed  in  a  temperature  controlled  environmental 
chamber  which  enabled  temperatures  to  be  maintained  within  ±1  °C.  The 
primary  component  and  user  interface  was  the  Hewlett-Packard  (HP) 
3562A  Dynamic  Signal  Analyzer  (DSA).  The  HP-3562A  was  used  to  provide 
a  swept  sine  signal  to  a  vibration  generator,  and  analyzed  the  returning 
data  signals.  The  HP-3562A  was  used  to  compute  the  frequency  response 
and  coherence  over  a  range  of  50  Hz  to  1050  Hz  using  the  discrete  Fourier 
Transform  in  the  swept  sine  mode.  Ten  averages  were  performed  at  each 
data  point  using  a  frequency  resolution  of  625  mHz  per  step.  The  source 
level  output  to  the  vibration  generator  was  set  at  1 .5  volts. 

A  schematic  of  the  testing  apparatus  is  shown  in  Figure  4.2.  Swept 
sine  source  signals  were  fed  from  the  output  jack  of  the  HP-3562A  to  a 
Wilcoxon  F3  vibration  generator  via  the  piezoelectric  output  of  a  Wilcoxon 
PA7C  power  amplifier.  The  vibration  generator  was  mounted  73.48  cm 
(28.93  in)  from  one  end,  and  12.7  cm  (5.0  in)  from  the  front  edge  of  each 
specimen  as  shown  in  Figure  4.3.  An  integral  force  transducer  was 
mounted  in  the  base  of  the  vibration  generator  to  measure  the  force  input  to 
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the  plate.  This  force  signal  was  then  fed  to  input  channel  one  of  the  DSA 
via  a  PCB  462-A  charge  amplifier.  Plate  accelerations  were  recorded  at 
various  points  using  a  PCB  303A-03  accelerometer  as  shown  in  Figure  4.3. 
Acceleration  data  was  fed  to  input  channel  two  of  the  DSA  via  a  PCB  482A05 
power  supply.  Frequency  response  and  coherence  data  was  then  recorded 
on  disk  for  further  analysis. 

Temperatures  within  the  testing  chamber  were  maintained  using  a 
NESLAB  RTE-8  refrigerated  circulating  bath  which  pumped  fluid  through 
a  small  heat  exchanger  in  the  testing  chamber  as  shown  in  Figure  4.2.  In 
order  to  accurately  monitor  the  temperature  of  the  plates,  a  small 
thermocouple  was  inserted  in  the  base  of  each  plate. 


SUPPORTING  BEAM 


Figure  4.1.  Configuration  of  plate  in  testing  chamber. 
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Figure  4.2.  Experimental  testing  system  schematic. 
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Figure  4.3.  View  of  shaker  and  accelerometer  locations. 
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B.  TESTING  PROCEDURE 


1 .  Undamped  reference  plate 

An  undamped,  reference,  frequency  response  measurement 
was  made  at  a  temperature  of  15.6  °C  (60  °F)  to  set  a  standard  response  by 
which  to  measure  the  effectiveness  of  the  damping  treatments.  The 
undamped  frequency  response  was  recorded  over  a  frequency  range  of 
50  - 1050  Hz  using  a  resolution  of  625  mHz  per  point  in  the  DSA. 

2.  Damped  plate  measurements 

Frequency  response  measurements  of  the  damped  plates  were 
made  at  a  temperature  of  15.6  °C  (60  °F)  at  several  nodes  on  the  plates  in 
order  to  capture  the  damped  response  of  as  many  modes  as  possible.  A 
representation  of  these  nodes  is  shown  in  Figure  4.3.  Responses  were 
recorded  over  a  frequency  range  of  50  - 1050  Hz  with  a  resolution  of 
625  mHz  per  point  in  the  DSA.  Zoom  measurements  were  also  made  to 
capture  better  data  for  certain  modes.  Modal  loss  factors  were  then 
estimated  from  the  frequency  response  and  coherence  measurements 
using  a  curve-fitting  technique  described  in  Reference  [14]. 

C.  SINGLE  DAMPING  LAYER  RESULTS 

The  single  damping  layer  treatment  was  tested  at  4.44  °C  (40  °F), 

15.6  °C  (60  °F),  and  26.7  °C  (80  °F)  so  that  the  effects  of  temperature  on  the 
damping  treatment  could  be  determined.  A  plot  of  the  single  damping 
layer  frequency  response  at  15.6  °C  is  shown  in  Figure  4.4.  The  single  layer 
damping  treatment  resulted  in  high  damping  with  modal  loss  factors 
ranging  from  0.223  at  53.4  Hz  to  0.091  at  876.6  Hz.  Due  to  the  coupling  of 
modes,  loss  factors  for  all  modes  were  not  measured.  The  frequency 
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response  of  the  damped  plate  is  characterized  by  a  frequency  shift  to  the  left 
and  a  smoothing  of  the  frequency  response  when  compared  to  the 
undamped  reference  plate.  The  single  layer  treatment  was  especially 
effective  at  reducing  the  frequency  response  of  a  mode  cluster  between  650 
and  950  Hz.  The  frequency  band  of  this  cluster  was  shifted  approximately 
200  Hz  with  the  amplitudes  of  the  responses  of  the  modes  being 
dramatically  reduced.  The  single  layer  treatment  was  also  effective  at 
reducing  the  amplitude  of  the  response  peaks  over  the  entire  spectrum  of 
measurement.  On  average,  the  highest  peaks  of  the  frequency  response  in 
the  undamped  condition  were  reduced  by  25  decibels,  a  reduction  of  17.8 
times. 

The  effect  of  temperature  on  the  damping  was  quite  pronounced  as 
shown  in  Figure  4.5.  As  the  testing  temperature  was  decreased,  the 
viscoelastic  layer  became  stiffer  and  damping  levels  were  increased.  A 
comparative  listing  of  the  loss  factors  at  different  temperature  is  in  Table 
4.1  and  a  plot  of  the  modal  loss  factors  is  shown  in  Figure  4.6.  Figure  4.5 
shows  the  trend  of  increased  damping  with  temperature  decrease,  and  a 
corresponding  shift  of  modal  frequencies  to  the  right  as  the  viscoelastic 
becomes  stiffer.  These  changes  are  especially  discernible  at  the  lower 
frequencies. 


Fxd  y  50  Hz  SINGLE  15.6  C  1.05k 


Figure  4.4.  Frequency  response  of  the  single  damping  layer  at  16.6  °C. 
t - :  undamped  response;  - :  damped  response] 


Figure  4.5.  Frequency  response  of  the  single  damping  layer  configuration 

at  different  temperatures. 

[ - :  4.44  °C, - :  15.6  °C, - :  26.7  °C] 
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TABLE  4.L  MODAL  LOSS  FACTORS  FOR  THE  SINGLE  DAMPING 
LAYER  AT  DIFFERENT  TEMPERATURE. 


4.44  °C  (40  °F)  15.6  °C  (60  °F)  26.7  °C  (80  °F) 


££Hzl 

a 

UHzl 

a 

£IHz) 

a 

64.3 

0.217 

53.7 

0.223 

49.0 

0117 

102.8 

0.223 

89.1 

0145 

83.5 

0.089 

157.4 

0.203 

138.9 

0172 

129.6 

0.082 

208.4 

0.183 

188.1 

0184 

177.9 

0.072 

327.5 

0.158 

301.6 

0120 

243.9 

0.062 

440.1 

0.144 

424.8 

0109 

409.4 

0.053 

474.9 

0.205 

476.8 

0111 

450.3 

0.066 

653.1 

0.079 

608.1 

0.067 

552.9 

0.068 

768.0 

0135 

644.2 

0.069 

586.7 

0.052 

852.0 

0.093 

722.5 

0130 

629.3 

0.042 

919.0 

0104 

817.6 

0.081 

681.4 

0.096 

953.5 

0.073 

876.6 

0.091 

790.1 

0.049 

858.2  0.056 

978.1  0.069 


0  200  400  600  800  1000 


Frequency  (Hz) 

Figure  4.6.  Comparison  of  modal  loss  factors  for  the  single  layer 
configuration  at  different  temperatures. 
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D.  DOUBLE  DAMPING  LAYER  RESULTS 


The  double  layer  damping  configuration  was  also  tested  at  4.44  °C, 
15.6  °C,  and  26.7  °C.  The  frequency  response  of  this  configuration  at  15.6  °C 
as  compared  to  the  undamped  reference  plate  is  shown  in  Figure  4.7. 
Damping  in  the  double  layer  configuration  is  also  high,  with  modal  loss 
factors  ranging  form  0.301  at  53.3  Hz  to  0.107  at  832.4  Hz.  Due  to  modal 
coupling  loss  factors  for  all  modes  were  not  measured.  The  frequency 
response  of  the  two-layer  configuration  is  also  characterized  by  a  dramatic 
reduction  in  response  amplitude  and  a  frequency  shift  to  the  left.  The  peak 
undamped  responses  were  reduced  by  an  average  of  27  decibels,  or  a 
reduction  of  22.4  times  from  the  reference  condition. 

The  effect  of  temperature  on  the  double  layer  damping  treatment  is 
shown  in  Figure  4.8.  As  with  the  single  layer  case,  damping  in  the  double 
layer  configuration  increased  with  a  decrease  in  temperature.  This 
configuration  also  shows  the  shift  of  modal  frequencies  to  the  right  as 
temperature  decreases  and  the  viscoelastic  becomes  stiffer.  Modal  loss 
factors  for  the  double  layer  configuration  are  listed  in  Table  4.2,  and  are 
plotted  for  comparison  in  Figure  4.9. 

To  compare  the  effectiveness  of  the  single  layer  and  double  layer 
configurations  their  frequency  responses  at  15.6  °C  are  plotted  in  Figure 
4.10  with  a  plot  comparing  modal  loss  factors  in  Figure  4.11.  The  two 
responses  are  quite  similar  although  the  double  layer  configuration  does 
show  an  increase  of  approximately  22  percent  in  modal  loss  factor. 
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FxdXY  50  Hz  DOUBLE  LAYER  15.6  C  1 . 05k 


Figure  4.7.  Frequency  response  of  the  double  layer  at  15.6  °C. 
[ - :  undamped  response  , - :  damped  response] 


FxdXY  50  Hz  4.44  /  15.6  /  26.7  C  1 . 05k 

Figure  4.8.  Frequency  response  comparison  of  the  double  layer 


configuration  at  different  temperatures. 

[ - :  4.44  °C, - :  15.6  °C  , - :  26.7  °C] 
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TABLE  4JL  MODAL  LOSS  FACTORS  FOR  THE  DOUBLE  LAYER 
CONFIGURATION  AT  DIFFERENT  TEMPERATURE 


4.44  °C  (40  °F) 


am 

n 

51.9 

0.273 

122.5 

0.224 

212.7 

0.187 

290.5 

0.197 

382.7 

0.174 

494.2 

0.144 

558.2 

0.169 

602.2 

0.198 

642.4 

0.157 

752.8 

0168 

815.4 

0159 

896.3 

0156 

968.3 

0130 

15.6  °C  (60  °F) 


£IHz) 

21 

55.3 

0.301 

871 

0.215 

142.1 

0.217 

190.5 

0.212 

297.8 

0.139 

366.8 

0.154 

419.3 

0.125 

441.2 

0100 

618.7 

0.098 

680.9 

0.096 

715.5 

0.060 

832.4 

0107 

26.7  °C  (80  °F) 


f  (Hz) 

n 

49.7 

0.202 

81.9 

0.188 

127.6 

0156 

172.4 

0.117 

278.2 

0.097 

317.7 

0.077 

396.1 

0.070 

428.5 

0.077 

602.0 

0.067 

638.4 

0.072 

738.0 

0.071 

802.6 

0.077 

8461 

0.050 

Figure  4.9.  Comparison  of  modal  loss  factors  for  the  double  layer  at 

different  temperature. 
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E.  POCKET  PLATE  RESULTS 


The  milled  pocket  plate  was  constructed  as  previously  shown  in 
Figure  3.3.  The  constraining  layer,  or  cover  plate,  was  welded  in  place 
using  tack  welds  in  an  attempt  to  keep  the  damping  material  away  from  the 
heat  of  welding,  and  the  cover  plate  from  warping,  instead  of  using  a 
continuous  weld  bead  as  was  previously  attempted  [Ref.  2].  The  cover  plate 
was  welded  to  the  base  at  the  comers,  at  the  midpoint  of  the  short  side  and 
at  three  equally  spaced  locations  along  the  long  dimension  as  shown  in 
Figure  4.12.  Following  welding  the  plate  was  tested  to  ensure  that  the 
viscoelastic  had  not  been  damaged  by  the  heat  of  welding. 


Figure  4.12.  Location  of  tack  welds  on  the  cover  plate 
of  the  pocket  plate  configuration. 

The  pocket  plate  was  tested  at  15.6  °C  (60  °F)  and  the  frequency 
response  is  shown  in  Figure  4.13.  The  response  indicates  that  the 
viscoelastic  layer  was  not  damaged  by  welding  and  that  good  damping  was 
attained.  Modal  loss  factors  ranged  firom  0.067  at  62.1  Hz  to  0.090  at  923  Hz. 
Although  damping  is  good,  it  is  approximately  half  that  of  the  single  layer 
configuration.  One  reason  for  this  is  that  the  viscoelastic  material  does  not 


completely  cover  the  base  structure.  Another  reason  is  the  presence  of  the 
welded  cover  plate.  Due  to  the  welded  conditions  the  cover  plate  cannot 
induce  shear  deformation  in  the  viscoelastic  layer  as  well  as  a  true 
constraining  layer,  and  therefore  produces  less  damping  than  the  single 
layer  configuration.  The  effects  of  the  welded  cover  plate  are  especially  felt 
in  modes  below  300  Hz  where  the  frequency  response  is  quite  peaked.  The 
response  curve  becomes  more  rounded  and  the  effects  of  the  damping  layer 
are  seen  as  frequency  increases. 

The  modal  loss  factors  for  the  pocket  plate  are  listed  in  Table  4.3  and 
are  plotted  in  Figure  4.14.  Even  though  the  damping  is  less  than  the  single 
layer,  the  plate  is  still  adequately  damped  as  shown  in  the  frequency 
response  plot  in  Figure  4.13.  In  this  configuration  the  modal  loss  factors 
remained  relatively  constant  throughout  the  testing  spectrum.  The 
increase  in  modal  loss  factor  values  above  800  Hz  is  due  primarily  to  modal 
coupling,  and  the  measured  modal  loss  factors  in  this  range  are  not 
reliable. 


Fxd  V  50  Hz  POCKET  PLATE  15.6  C  1 . 05k 


Figure  4.13.  Frequency  response  of  the  pocket  plate  at  15.6  °C  (60  °F). 
[ - :  reference  plate  , - :  pocket  plate  ] 
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Modal  Loss  Factors 


TABLE  43  MODAL  LOSS  FACTORS  FOR  THE 
POCKET  PLATE  CONFIGURATION. 


f  (Hz) 

21 

62.1 

0.067 

93.4 

0.042 

142.0 

0.056 

194.9 

0.060 

308.9 

0.051 

437.5 

0.041 

495.6 

0.036 

HHz) 

n 

565.3 

0.043 

625.7 

0.047 

643.2 

0.049 

687.0 

0.044 

841.9 

0.056 

891.5 

0.077 

923.0 

0.090 

Figure  4.14.  Modal  loss  factors  versus  frequency  for 
the  pocket  plate  configuration. 
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FLOATING  ELEMENT  RESULTS 


The  milled  “floating  element”  configuration  was  constructed  as 
previously  shown  in  Figure  3.7.  The  cover  plate  was  welded  in  a  fashion 
similar  to  the  pocket  plate  as  shown  in  Figure  4.12.  The  center 
constraining  layer,  or  floating  element,  was  made  slightly  smaller  than  the 
surrounding  structure  thus  allowing  the  floating  element  to  act  as  a  “true” 
constraining  layer. 

The  frequency  response  of  the  floating  element  configuration  at 
15.6  °C  (60  °F)  is  shown  in  Figure  4.15.  The  floating  element  is  quite 
effective  as  the  response  shows  a  good  reduction  in  peak  modal  response. 
Measured  modal  loss  factors  range  from  0.089  at  66  Hz  to  0.064  at  935  Hz.  A 
listing  of  measured  modal  loss  facto  rs  is  in  Table  4.4  and  are  plotted  in 
Figure  4.16.  As  with  the  previous  cases,  the  frequency  response  of  the 
floating  element  configuration  is  characterized  by  a  frequency  shift  to  the 
left  and  a  smoothing  of  the  response  as  frequency  increases. 

In  a  comparison  of  the  pocket  plate  and  floating  element 
configurations,  the  frequency  responses  are  plotted  in  Figure  4.17.  A 
comparison  of  modal  loss  factors  for  the  ‘wo  configurations  is  shown  in 
Figure  4.18.  The  two  frequency  response  plots  are  similar,  however,  the 
frequency  response  of  the  floating  element  configuration  is  more  rounded 
than  that  of  the  pocket  plate.  The  major  difference  between  the  two 
configurations  is  seen  in  Figure  4.18.  Modal  loss  factors  for  the  floating 
element  show  an  average  increase  of  25  percent  over  those  of  the  pocket 
plate.  Reasons  for  this  increase  are  the  added  constraining  effect  of  the 
floating  element  and  additional  layer  of  damping  material  present 
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Modal  Loss  Factors 


TABLE  4.4  MEASURED  MODAL  LOSS  FACTORS  FOR  THE  FLOATING 

ELEMENT  CONFIGURATION. 


mi zi 

n 

f  (Hz) 

n 

66.0 

0.089 

508.1 

0.089 

104.9 

0.040 

540.6 

0.115 

154.0 

0.094 

618.7 

0.121 

211.4 

0.058 

656.7 

0.100 

272.8 

0.063 

708.6 

0.088 

323.5 

0.075 

857.6 

0.089 

435.8 

0.090 

935.4 

0.064 

0.14 
0.12 
0.10 
0.03 
0.06 
0.04 
0.02 

0  200  400  600  800  1000 

Frequency  (Hz) 

Figure  4.16.  Modal  loss  factors  for  the  floating  element  configuration. 
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FREQ  RESP 


Figure  4.17.  Frequency  response  comparison  for  the  pocket  plate  and 
floating  element  configurations. 

[ - :  pocket  plate  , - :  floating  element  ] 


Frequency  (Hz) 

Figure  4.18.  Modal  loss  factor  comparison  for  the  pocket  plate  and  floating 

element  configurations. 
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V.  FINITE  ELEMENT  RESULTS 


A.  UNDAMPED  REFERENCE  PLATE 

The  first  step  in  the  finite  element  analysis  procedure  was  to  model 
and  analyze  the  undamped  reference  plate  for  its  modal  frequencies  and 
frequency  response.  The  finite  element  model  was  generated  using 
PATRAN,  a  computer  aided  interactive  graphics  program  developed  by 
PDA  Engineering.  PATRAN  is  widely  used  for  conceptual,  preliminary, 
and  detailed  design  and  analysis  of  complex  systems.  One  of  PATRAN’s 
major  advantages  is  in  the  interactive  construction  of  finite  element  models 
for  use  by  MSC/NASTRAN,  and  its  ability  to  display  MSC/NASTRAN 
results  in  an  easily  understood  graphic  format  [Ref.  15]. 

The  reference  plate  was  modeled  using  84  plate  (QUAD4)  elements  as 
shown  in  Figure  5.1.  The  QUAD4  element  is  an  isoparametric  element 
with  four  nodes,  one  at  each  comer  of  the  element  [Ref.  7].  A  normal  mode 
extraction  was  then  performed  in  order  to  compare  numerical  results  with 
experimentally  determined  modal  frequencies.  Once  satisfactory 
agreement  between  the  modal  frequencies  calculated  in  NASTRAN  and 
those  obtained  experimentally  was  attained,  a  direct  frequency  response 
calculation  was  performed  in  NASTRAN.  This  frequency  response  was 
then  used  as  the  reference  response  for  further  finite  element  models 
incorporating  viscoelastic  damping  treatments.  The  excitation  for  the 
frequency  response  calculation  was  a  sinusoidal  force  with  an  amplitude  of 
1 .0  applied  at  node  66.  The  response  point  was  node  58  as  shown  in  Figure 
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5.1.  These  two  nodes  correspond  to  the  points  on  the  plates  used  in  the 
experimental  portion  of  the  research  where  the  vibration  generator  and 
accelerometer  were  attached. 


Figure  5.1.  Finite  element  model  of  the  undamped  reference  plate. 

B.  SINGLE  DAMPING  LAYER 

The  modeling  of  the  single  constrained  layer  damping  system  was 
done  using  techniques  described  by  Job  *on  and  Kienholz  [Ref.  1].  As 
shown  in  Figure  5.2,  the  viscoelastic  layer  is  modeled  using  solid  (HEXA) 
elements,  while  the  base  layer  and  constraining  layer  were  modeled  using 
QUAD4  elements.  The  HEXA  element  is  a  solid,  isoparametric  element 
having  eight  nodes,  one  at  each  corner  of  the  element  with  three 
translational  degrees  of  freedom  at  each  node  [Ref.  7].  The  use  of  solid 
elements  for  the  viscoelastic  layer  allows  the  strain  energy  due  to  shearing 
to  be  adequately  represented.  Plate  elements  are  used  in  the  base  layer  and 
constraining  layer  because  of  their  ability  to  account  for  stretching  and 


bending  deformations.  The  plate  element  allows  its  nodes  to  be  offset  from 
the  plate’s  center  to  the  surface  of  the  plate,  coincident  with  the  corner 
nodes  of  the  solid  viscoelastic  elements  [Ref.  1].  Thus,  the  single 
constrained  layer  system  was  modeled  using  only  two  layers  of  nodes,  a 
simple  process  in  PATRAN.  For  the  single  layer  configuration,  a  model 
having  84  elements  per  layer,  and  an  element  meshing  scheme  shown  in 
Figure  5.2  was  used. 


Figure  5.2.  Finite  element  representation  of  the 
single  constrained  layer. 


Once  the  damped  plate  had  been  modeled,  normal  mode  extractions 
were  made  using  MSC/NASTRAN.  Five  separate  runs  were  conducted 
using  the  material  properties  of  ISD-112  at  50,  200,  500,  800,  and  1000  Hz 
and  at  a  temperature  of  15.6  °C  (60  °F).  In  addition  to  the  modal 
frequencies,  the  strain  energy  in  the  viscoelastic  elements  and  the  entire 
model  were  output  from  NASTRAN. 

Since  the  shear  modulus  of  a  viscoelastic  material  changes  with 
frequency,  it  was  necessary  to  estimate  the  actual  modal  frequencies  of  the 
damped  plate  using  an  interpolation  procedure  outlined  by  Johnson  and 
Kienholz  [Ref.  16].  The  first  step  of  the  interpolation  process  was  to  plot  the 
shear  modulus  of  ISD-112  versus  frequency  from  5  to  1000  Hz.  Then,  for 
the  first  mode,  using  NASTRAN  results  based  on  ISD-112  material 
properties  a  50  Hz,  the  first  modal  frequency  predicted  by  NASTRAN  and 
the  corresponding  shear  modulus  were  plotted.  The  same  was  then  done 
using  the  first  natural  frequency  predicted  by  normal  mode  extraction 
based  on  ISD-112  material  properties  at  frequencies  of  500  and  1000  Hz.  A 
curve  was  then  passed  through  these  three  points.  The  point  where  the 
NASTRAN  modal  frequencies  for  the  first  mode  intersected  the  ISD-112 
shear  modulus  curve  was  taken  to  be  the  interpolated  modal  frequency  of 
the  single  damping  layer  configuration.  A  plot  of  the  intersection  of  these 
two  curves  for  the  first  and  second  modes  is  shown  in  Figure  5.3.  This 
interpolation  process  was  then  repeated  for  each  mode  through  1000  Hz. 


Figure  5.3.  Interpolation  of  the  first  and  second  modal  frequencies 
for  the  single  constraining  layer  configuration. 

Once  the  interpolated  modal  frequencies  were  found,  the  modal 
strain  energy  equations  (2.14)  and  (2.15)  were  used  to  compute  modal  loss 
factors  for  the  single  layer  configuration.  A  set  of  modal  loss  factors  was 
computed  based  on  the  modal  strain  energies  computed  using  viscoelastic 
properties  at  reference  frequencies  of  50  ,200,  500,  and  800  Hz.  A  set  of 
composite  modal  loss  factors  for  the  modal  frequencies  near  these  reference 
frequencies  was  then  selected.  The  resulting  modal  loss  factors  are  shown 
in  Table  5.1  and  are  plotted  versus  frequency  in  Figure  5.4.  As  seen  in 
Figure  5.4,  the  modal  strain  energy  method  is  predicting  high  damping  for 
this  configuration  with  an  average  modal  loss  factor  of  0.195. 
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Modal  Loss  Factor 


TABLE  5.L  MODAL  LOSS  FACTORS  FOR  THE  SINGLE 
LAYER  CONFIGURATION  AS  COMPUTED  USING 
THE  MODAL  STRAIN  ENERGY  METHOD 


f  (Hz) 

21 

f  (Hz) 

n 

63 

0.249 

517 

0.237 

95 

0.181 

586 

0.217 

153 

0.234 

632 

0.164 

202 

0.212 

663 

0.155 

285 

0.229 

712 

0.180 

322 

0.214 

827 

0.134 

447 

0.200 

869 

0.152 

463 

0.108 

883 

0.115 

483 

0.258 

Figure  5.4.  Estimated  modal  loss  factors  for  the  single  layer  configuration 
with  the  curve  fit  used  for  the  MSC/NASTRAN  damping  table. 


Using  the  set  of  composite  modal  loss  factors,  the  modal  frequency 
response  of  the  damped  plate  was  computed  using  MSC/NASTRAN.  Modal 
damping  was  introduced  to  the  model  using  the  SDAMP  option  in  the  Case 
Control  Deck  and  the  TABDMP1  damping  table  in  the  Bulk  Data  Deck  as 
described  in  the  MSC/NASTRAN  Handbook  for  Dynamic  Analysis  [Ref.  17]. 
Since  NASTRAN  uses  a  linear  interpolation  between  points  in  the  damping 
table  to  describe  the  modal  damping  in  the  model  [Ref.  7],  a  simple  curve  fit 
was  applied  to  the  set  of  composite  modal  loss  factors  as  shown  in 
Figure  5.4.  Points  from  this  curve  fit  were  then  used  in  the  NASTRAN 
damping  table.  To  compute  the  modal  frequency  response,  a  unit  excitation 
force  was  applied  at  the  same  node  as  the  undamped  plate,  and  the  node 
used  for  the  response  was  also  the  same  as  the  undamped  plate. 

The  results  of  the  modal  frequency  response  calculations  are  shown 
in  Figure  5.5.  The  dashed  line  represents  the  undamped  reference  plate, 
and  the  solid  line  represents  the  modal  frequency  response  of  the  single 
layer  configuration.  Material  properties  at  200  Hz  were  used  for  the 
ISD-112  damping  material.  The  first  thirty  modes  were  used  in  the  modal 
summation  for  the  response.  A  listing  of  the  MSC/NASTRAN  data  deck 
used  to  compute  the  modal  frequency  response  is  in  Appendix  D. 
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Figure  5.5.  Modal  frequency  response  of  the  single  layer 
configuration  as  computed  in  NASTRAN. 

[ - :  reference  plate  , - :  single  layer] 
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The  modal  loss  factors  estimated  using  the  modal  strain  energy 
method  are  compared  to  those  measured  experimentally  for  the  single  layer 
configuration  in  Figure  5.6.  The  estimated  loss  factors  are  greater  than  the 
experimentally  determined  loss  factors  throughout  the  spectrum  of 
interest,  and  especially  in  the  middle  frequencies. 

The  modal  frequency  response  of  the  single  layer  configuration  is 
compared  to  the  experimentally  measured  frequency  response  in 
Figure  5.7.  The  comparison  was  accomplished  by  normalizing  both  the 
experimentally  determined  frequency  response  and  the  frequency  response 
computed  in  NASTRAN.  Both  responses  were  normalized  using  a  value  of 
1.0  ^n^C'  •  The  effects  of  the  greater  loss  factors  estimated  by  the  modal 
strain  energy  method  are  obvious  as  the  level  of  the  predicted  response  is 
lower  than  the  measured  response.  The  shift  in  frequency  between  the  two 
curves  is  due  to  the  finite  element  model  being  inherently  stiffer  than  the 
actual  system.  The  correlation  between  the  two  curves  is  especially  good 
below  250  Hz  as  this  is  where  the  differences  between  estimated  and 
measured  modal  loss  factors  are  the  least. 
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Figure  5.6.  Comparison  of  estimated  and  measured  modal  loss 
factors  for  the  single  layer  configuration. 
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Figure  5.7.  Frequency  response  comparison  for 
the  single  layer  configuration. 
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C.  DOUBLE  DAMPING  LAYER 


The  modeling  of  the  double  constrained  layer  damping  system  was 
accomplished  as  shown  in  Figure  5.8.  The  double  layer  configuration 
consists  of  a  base  layer  modeled  with  offset  QUAD4  elements,  two 
viscoelastic  layers  consisting  of  HEXA  elements,  and  the  top  constraining 
layer  modeled  with  offset  QUAD4  elements.  The  middle  constraining  layer 
was  modeled  using  three  layers  of  HEXA  elements  in  order  to  give  this 
layer  the  stiffness  necessary  to  act  as  a  constraining  layer.  The  model  was 
meshed  using  60  elements  in  each  layer  as  shown  in  Figure  5.8. 

Using  this  model,  the  modal  strain  energy  method  was  applied  to 
determine  approximate  modal  frequencies  and  loss  factors  for  the  double 
layer  configuration.  To  determine  the  loss  factors,  normal  mode 
extractions  were  performed  using  reference  frequencies  of  50,  200,  500,  800, 
and  1000  Hz.  A  composite  set  of  modal  loss  factors  for  the  double  layer 
system  was  then  compiled  based  on  the  estimated  modal  frequency’s 
relation  to  the  reference  frequency  used  to  calculate  modal  strain  energies. 
This  composite  set  of  modal  loss  factors  is  listed  in  Table  5.2,  and  is  plotted 
in  Figure  5.19  as  a  comparison  to  the  experimentally  determined  modal 
loss  factors  for  the  double  layer  configuration.  The  estimated  modal  loss 
factors  for  the  double  layer  show  high  damping,  but  they  compare  favorably 
with  those  measured  experimentally. 
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TABLE  5.2.  ESTIMATED  MODAL  LOSS  FACTORS  FOR  THE  DOUBLE 

LAYER  CONFIGURATION. 
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Figure  5.9.  Modal  loss  factors  for  the  double  layer  configuration 
as  determined  from  the  modal  strain  energy  method 
and  determined  experimentally. 
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The  modal  frequency  response  of  the  double  layer  configuration  was 
computed  in  a  manner  similar  to  the  single  layer  in  that  smoothed  loss 
factor  data  was  used  in  the  MSC/NASTRAN  damping  table.  Likewise,  a 
unit  excitation  force  was  applied,  and  the  first  thirty  modes  were  used  in 
the  modal  summation.  The  results  of  the  modal  frequency  response 
calculation  are  shown  in  Figure  5.10.  The  dashed  line  represents  the 
undamped  response,  and  the  solid  line  represents  the  frequency  response  of 
the  double  layer  configuration. 

The  frequency  response  calculated  using  NASTRAN  was  compared 
to  the  experimentally  determined  frequency  response  of  the  double  layer 
configuration  as  shown  in  Figure  5.11.  Comparison  of  the  two  frequency 
response  curves  shows  similarity  in  form,  but  a  much  lower  response  level 
for  the  numerically  determined  response.  Once  again  this  could  be  due  to 
the  higher  damping  predicted  by  the  modal  strain  energy  method  and  the 
inherently  higher  stiffness  of  the  finite  element  model. 
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Figure  5.10.  Modal  frequency  response  for  the  double  layer 
configuration  as  computed  in  NASTRAN. 

[ - :  reference  plate  , - :  double  layer] 
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Figure  5.11.  Comparison  of  the  experimentally  determined 
and  numerically  predicted  frequency  responses 
for  the  double  layer  configuration. 
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D.  POCKET  PLATE  RESULTS 


The  pocket  plate  configuration  was  modeled  with  the  same  offset 
plate  elements  and  solid  viscoelastic  elements  as  the  single  layer 
configuration.  However,  the  pocket  plate  required  the  modeling  of  the 
milled  structure  around  the  viscoelastic  material  and  the  welds  between 
the  cover  plate  and  milled  plate.  A  representation  of  the  model  is  shown  in 
Figure  5.12.  The  base  structure,  cover  plate,  and  the  structure  immediately 
around  the  cover  plate  was  modeled  using  offset  plate  elements.  The 
viscoelastic  material  and  the  portion  of  the  structure  immediately  adjacent 
to  it  were  modeled  using  the  solid  HEXA  elements.  Since  the  viscoelastic 
material  and  cover  plate  are  physically  separated  from  the  surrounding 
plate,  except  where  the  viscoelastic  is  adhered  to  the  base  structure,  care 
was  necessary  in  creating  the  finite  element  mesh. 

The  model  was  created  in  PATRAN  using  PATRAN’s  node  editing 
and  equivalencing  capabilities  [Ref.  18].  This  allowed  the  generation  of  a 
finite  element  mesh  with  two  nodes  at  the  same  geometric  point  in  space. 
Using  this  node  editing  technique,  a  mesh  was  created  which  allowed  the 
viscoelastic  and  cover  plate  to  vibrate  separately  from  the  surrounding 
structure,  yet  at  the  same  time,  keep  the  number  of  elements  and  nodes  in 
the  model  to  a  minimum.  The  welded  points  on  the  cover  plate  were  also 
modeled  using  node  editing  techniques.  At  weld  points  the  finite  element 
node  on  the  cover  plate  was  equivalenced  with  its  corresponding  node  on  the 
base  structure,  resulting  in  a  single  node  and  a  connection  between  an 
otherwise  separate  base  structure  and  cover  plate.  At  non-welded  points 
on  the  cover  plate  there  were  two  nodes  at  the  same  geometric  point;  one  to 


represent  the  cover  plate,  and  the  other  to  represent  the  base  structure.  The 
model  was  meshed  using  a  5x11  mesh  resulting  in  40  elements  in  the  cover 
plate,  viscoelastic  and  base  layer  as  shown  in  Figure  5.12. 

Using  this  modeling  scheme,  the  modal  strain  energy  method  was 
employed  to  estimate  the  modal  frequencies  and  loss  factors.  Normal  mode 
extractions  were  made  using  viscoelastic  material  properties  at  50,  200,  500, 
800,  and  1000  Hz.  Using  these  reference  frequencies  a  composite  set  of 
modal  loss  factors  was  obtained.  These  loss  factors  are  listed  in  Table  5.3 
and  are  plotted  versus  frequency  in  Figure  5.13.  The  estimated  loss  factors 
give  good  damping  over  the  spectrum  of  interest  with  an  average  modal  loss 
factor  of  0.075.  As  with  the  previous  cases,  damping  values  used  for  the 
modal  frequency  response  calculation  came  from  a  curve-fit  to  the  set  of 
composite  modal  loss  factors. 

The  modal  frequency  response  of  the  pocket  plate  was  computed 
using  the  first  30  modes  and  viscoelastic  material  properties  at  200  Hz.  The 
resulting  estimated  frequency  response  is  shown  in  Figure  5.14.  The 
response  shows  a  definite  frequency  shift  to  the  left  along  with  good 
damping  of  the  frequency  response  when  compared  to  the  undamped 
response. 
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TABLE  5.3.  ESTIMATED  MODAL  LOSS  FACTORS  FOR  THE  POCKET 


PLATE  CONFIGURATION. 
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Figure  5.13.  Estimated  modal  loss  factors  versus 
frequency  for  the  pocket  plate. 
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for  the  pocket  plate  configuration. 
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The  estimated  modal  loss  factors  and  modal  frequency  response  for 
the  pocket  plate  were  compared  to  those  measured  experimentally.  The 
loss  factor  comparison  is  shown  is  Figure  5.15  and  the  frequency  response 
comparison  is  shown  in  Figure  5.16.  The  estimated  modal  loss  factors  are 
higher  than  those  measured  experimentally,  however,  the  frequency 
responses  compare  quite  favorably  with  each  other.  The  frequency 
response  curve  calculated  through  finite  element  analysis  has  a  lower 
response  level  and  a  frequency  shift  to  the  right  of  the  measured  frequency 
response.  This  is  expected  due  to  the  increase  in  damping  predicted  by  the 
modal  strain  energy  method  and  by  the  fact  that  the  finite  element  model  is 
inherently  stiffer  than  the  physical  system. 
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Figure  5.15  Comparison  of  experimentally  measured  and  numerically 
estimated  modal  loss  factors  for  the  pocket  plate  configuration. 
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Figure  5.16.  Comparison  of  experimental  and  predicted  frequency 
responses  for  the  pocket  plate  configuration. 
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VL  CONCLUSIONS 


Constrained  viscoelastic  layer  damping  is  an  extremely  effective 
method  for  reducing  broadband  vibration.  In  each  of  the  experimental 
cases  the  peak  amplitudes  of  frequency  response  were  reduced  by 
approximately  25  decibels,  a  reduction  of  18  times  below  the  undamped 
reference  plate. 

Of  particular  note  are  the  performance  of  the  pocket  plate  and 
floating  element  configurations.  Although  they  are  not  ideal 
configurations  in  terms  of  “true”  constrained  viscoelastic  layer  damping, 
the  damping  levels  achieved  are  quite  satisfactory.  In  a  final  comparison  of 
experimentally  determined  modal  loss  factors  for  the  four  treatments,  the 
double  layer  configuration  yields  the  greatest  damping  and  the  pocket  plate 
the  least.  This  comparison  is  shown  in  Figure  6.1.  As  was  previously 
reported  in  Section  IV,  tl  floating  element  configuration  yielded  an 
average  increase  of  25  percent  over  the  modal  loss  factors  of  the  pocket 
plate.  It  is  also  noted  that  the  average  modal  loss  factor  for  the  pocket  plate 
is  approximately  50  percent  of  the  average  modal  loss  factor  of  the  single 
layer  configuration.  Similarly,  the  average  modal  loss  factor  for  the 
floating  element  configuration  is  approximately  50  percent  that  of  the 
double  layer  treatment. 
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Figure  6.1.  Experimental  modal  loss  factor  comparison 
for  the  four  damping  configurations. 


The  modal  strain  energy  method  tends  to  overpredict  the  modal  loss 
factors  for  these  highly  damped,  thick  plates.  Although  the  modal  strain 
energy  method  predicted  modal  loss  factors  greater  than  those  measured, 
the  same  relative  differences  in  modal  loss  factor  between  damping 
configurations  are  maintained  as  shown  in  Figure  6.2.  The  estimated 
modal  loss  factors  for  the  pocket  plate  are  approximately  60  percent  less 
than  those  of  the  single  layer  configuration  This  indicates  that  there  is  a 
consistency  between  the  damping  values  predicted  by  the  modal  strain 
energy  method  and  those  of  the  actual  physical  system. 
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Figure  6.2.  Comparison  of  numerically  estimated  modal  loss  factors  for 
the  single  layer,  double  layer,  and  pocket  plate  configurations. 
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There  are  several  possible  reasons  for  the  differences  between  the 
experimentally  determined  modal  loss  factors  and  those  estimated  by  the 
modal  strain  energy  method.  The  first  is  that  the  material  properties  of 
ISD-112  as  reported  by  the  3M  Corporation  on  the  reduced  frequency 
nomogram  may  not  be  consistent  with  the  material  properties  actually 
present  in  the  material  used.  Since  the  numerical  analysis  was  based  on 
the  reported  material  properties  this  is  a  possible  source  of  uncertainty  in 
the  results. 

The  second  source  of  uncertainty  is  in  the  adhesion  of  the  ISD-112  to 
the  aluminum  plates.  Although  the  plates  were  clean  when  the  ISD-112 
was  applied,  it  was  noted  that  the  adhesive  qualities  of  the  ISD-112  were  not 
uniform  throughout  the  material.  A  lack  of  adhesion  may  cause  a 
reduction  in  the  damping  capability  of  the  system. 


Another  source  for  the  difference  between  the  experimental  results 
and  numerical  results  lies  within  the  finite  element  model.  The  number  of 
elements  used  in  the  model  greatly  affects  its  “stiffness.”  As  the  number  of 
elements  is  increased  the  model  should  become  less  stiff  and  results  are 
expected  to  approach  those  that  are  measured.  Also,  the  type  of  element 
used  to  model  the  base  layer  and  constraining  layer  may  have  an  effect.  In 
this  research  plate  elements  were  used  to  model  both  the  base  and 
constraining  layers.  It  is  possible  that  one  or  more  layers  of  solid  elements 
may  produce  results  that  agree  better  with  experimental  results. 

The  modal  strain  energy  method  is  a  useful  tool  for  vibration 
damping  design,  however,  the  method  appears  to  overestimate  modal  loss 
factors  for  thick,  highly  damped  plates  by  as  much  as  50  percent.  In  each 
configuration  that  was  modeled,  the  modal  loss  factors  estimated  below 
300  Hz  were  fairly  close  to  the  actual  values,  and  above  300  Hz  the  difference 
between  estimated  as  measured  loss  factors  increased. 

One  drawback  to  the  modal  strain  energy  method  in  design  is  the 
CPU  time  required  for  normal  mode  extraction  and  modal  frequency 
response;  especially  in  complex  structures  with  a  large  number  of 
elements.  Although  the  modal  strain  energy  method  is  good  for  finalizing 
a  design,  it  should  not  be  used  for  design  optimization  due  to  the  amount  of 
CPU  time  required  for  normal  mode  extraction  and  modal  frequency 
response  calculations. 
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VIL  RECOMMENDATIONS 


The  thick  plate  used  in  this  research  is  a  generic  model  of  many 
physical  systems  that  may  see  use  in  naval  application.  There  are  several 
areas  which  deserve  more  research  and  clarification,  including  the 
following: 

•  Model  the  floating  element  configuration  in  finite  elements  and 
check  the  effectiveness  of  the  modal  strain  energy  method  in 
predicting  modal  loss  factors  for  this  configuration. 

•  Investigate  the  capability  of  the  modal  strain  energy  method  to 
predict  loss  factors  for  thick,  lightly  damped  plate  configurations.  It 
is  possible  that  the  modal  strain  energy  method  works  better  for 
lightly  damped  configurations. 

•  Investigate  methods  for  improving  the  adhesion  of  ISD-112  to  the 
base  structure  and  constraining  layer. 

•  Investigate  the  possible  use  of  the  floating  element  as  an  additional 
vibration  absorber  by  changing  its  mass  in  order  to  damp  a  specific 
frequency.  In  this  same  idea,  the  use  of  tuned  dampers  along  with 
the  constrained  layer  treatments  could  be  looked  at. 


APPENDIX  A 


FORTRAN  PROGRAM  USED  TO  COMPUTE  MODAL  LOSS  FACTORS 
FOR  THE  SINGLE  DAMPING  LAYER  DESIGN 

This  program  used  the  Ross-Kerwin-Ungar  equations  of  Section  III 
to  compute  estimated  modal  frequencies  and  loss  factors  for  the  single  layer 
configuration.  Modal  frequencies  for  an  undamped  plate  are  read  from  a 
data  file  and  estimated  modal  frequencies  and  loss  factors  for  various  layer 
thicknesses  are  output  to  another  file.  Material  properties  of  ISD-112  are 
computed  using  University  of  Dayton  data  and  curve-fitting  equations  to  the 
reduced  frequency  nomogram  [Ref.  9,12].  The  units  used  in  this  program 
are  pounds,  inches,  and  seconds. 
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PROGRAM  CARPET 

XXXXXXXyyXXKXXXXXXXXXHXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
THIS  PROGRAM  IP,  TO  CALCULAI  E  SYSTEM  LOSS  FACTORS  POR  VARIOUS 
T  EHP  ER  A  I  URE  S  Alin  MOOES  OF  A  COIIS  I  RAIIIEI)  LAYER  VISCOELASTIC 
DAMPING  SYSTEM  USING  THE  ROSS- K ERH 1 II- UIIGAR  EQUATIONS.  SYSTEM  LOSS 
TACIORS  ARE  COMPUIPU  TOR  A  I  EMI'  E  R A I  UR  E  RANGE  OP  30-100  DEGREES 
FAREIIUE I T .  BASE  rL AT  E  THICKNESSES  VARY  FROM  0.375  TO  n.75  INCHES 
III  0.125  1 NCREMEII !  3  .  FOR  EACH  BASE  PLATE  THICKNESS,  THE  VISCOFLASTIC 
THICKNESS  IS  VARIED  TROM  0.015  TO  0.060  INCHES  111  0.015  INCH 
INCITEMENTS.  FOR  EACH  VISCOELASTIC  THICKNESS  THE  CONSTRAINING 
LAYER  THICKNESS  IS  VARIED  FROM  0.0625  TO  0.25  INCHES  III 
0  0625  INCH  INCREMENTS. 

THIS  PROGRAM  APLIES  TO  A  FREE-FREE-FREE-TREE  PLATE. 

THE  VISCOELASTIC  MATERIAL  IS  3M  ISD-112. 

THE  FULL  Oil  I  till  COFFFICI  F.HTS  ARC  DEFINED: 

El  =  YOUNG'S  MODULUS  OF  BASE  PLATE  (PSI) 

E 2  =  YOUNG'S  MODULUS  OF  VISCOELASTIC  LAYER  (PSI) 

C 3  =  YOUNG'S  MODULUS  OF  CONSTRAINING  LAYER  (PSI) 

02  -  SHEAR  MODULUS  OF  VISCOELASTIC  MATERIAL  (PSI) 

Hill  =  POISSON'S  RATIO  OF  BASE  PLATE 
IIU2  -  POISSON'S  RATIO  OF  VISCOELASTIC  MATERIAL 

111  -  THICKNESS  OF  BASE  PLATE  (IN) 

112  =  THICKNESS  OF  VISCOELASTIC  LAYER  (  1 N ) 

113  =  THICKNESS  OF  CONSTRAINING  LAYER  (IN) 

HINT  =  TOTAL  PLATE  THICKNESS  (IN) 

RII01  =  DENSITY  OF  BASE  PLATE  (LBF-SEC**2/IN**A) 

RII02  =  DENSITY  OF  VISC0ELAS1TC  MATERIAL 
RH03  =  DENSITY  OF  CONSTRAINING  LAYER 

I  0  ,  T  ROM ,  MROM  .  IT,  ML  ,  ET  FROL  ,  SL  ,  Stl ,  TROI.  ,  8  C  ARE  COEFFICIENTS 
FOR  THE  REDUCED  FREQUENCY  NOMOGRAM  EQUATIONS 
FP  -  MODAL  FREQUENCY  OF  THE  UNDAMPED  BASE  PLATE  (HZ) 

TCP  =  MODAL  FREQUENCY  OF  THE  COMPOSITE  PLATE  (HZ) 

E I  A  2  =  LOSS  FACTOR  OF  VISCOELASTIC  MATERIAL 
ETAS  =  SYSTEM  LOSS  FACTOR 
KQR  -  HAVE  NUMBER  OF  PLATE 
GC  =  GRAVITATIONAL  CONSTANT  (IN/SEC**2) 

I  =  TEMPERATURE  OT  VISCOELASTIC  MATERIAL  (DEG  F) 

HP  =  MODAL  FREQUENCY  OF  UNDAMPED  PLATE  (RAD/SEC) 

TP  -  MODAL  FREQUENCY  OF  UNDAMPED  PLATE  (HZ) 

HIT  =  MODAL  FREQUENCY  OF  DAMPED  PIAIE  (RAD/SEC) 

TCP  =  MODAL  FREQUENCY  OF  DAMPED  PLATE  (HZ) 

THE  UNITS  USED  IN  THIS  PROGRAM  ARE  LB,  INCH,  SEC,  DEGREES  F 

UNDAMPED  MODAL  FREQUENCIES  ARE  FROM  A  FINITE  ELEMENT  ANALYSIS 
OF  THE  TREE-FREE-FREE-FREE  FLATE . 

RESULTS  APE  OUTPUT  TO  DATA  FILE  "LOSFCTR  DATA" 

PRIOR  TO  PUNNING  THE  PROGRAM  TYPE  THE  COMMAND 
"FILEDEF  LOSTCTR  DISK  LOSTCTR  DATA" 

FINITE  El  EMC NT  MODAL  FREQUENCIES  ARE  INPUT  FROM  TILE: 

"PL ITRQ  DA  I  A" 

XXX vxx*X*XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

real  El , C2, E3, NUI , NU2, Ml , H2, H3, RNOI . RH02 , RH03 , HTOT 
REAL  10,  [  ROM, MROM, N, ML , ET T ROL , SI , SH, FROL , C 
REA  I  TP,  TCP, ETA2, ETAS, KQR, PI ,GC, T 
RE A  l  TRIO, ETA210, 5UB1, SUB2,SUB3,M10 
REAL  (12  1  ,  H31  ,  G.  Cl  ,  ALPHRE.  AirillM,  DRE,  DIM,  DEL  RE 
REAL  DEI  111,  EHCUBE ,  5UB6 ,  SUB5 ,  IIP , HCP ,  DENS ,  SUB6  ,  SUB7  ,  SUBS 
I  III  EGER  V 
DIMENSION  LIFTS) 

T I  r «* .  x A 1  AIK  1  .  ) 

OPEIK  IIIII  r-1  0  ,  ril  E-  '  L0SFC1  R  '  ,  ST  AT  US  -  '  01  D'  ) 
OPENCUNI!=ll,riLL=TLirRQ',SlAlUS='OLD') 

C  ASSIGN  MATERIAL  C0IISTAN1S  FOR  BASE  r  L  AT  E  ,  CONST  R  A I II I  TIG  LAYER 
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C  Ann  VISCOELASTIC 

c 

F  1  =  1  .  0E7 
HI =0 . 750 
RHO 1 =  .  0168 
Mill  =  .  33 
RHUc - . 035 
E  3  = 1  .  0E7 
RH03- .0968 
OC =386 . 

NU2  = . 5 
C 

C  DCrillE  C0IISTAIIT3  FOR  NOMOGRAM  EQUATIONS 

C 

T 0  =  109 . 0 
FROM =2 . 0F9 
HP.0M  =  688 .99 
11=  .  27 5 
ML  =3 . 7 
El  FROL  =1  .  08 
51=0.95 
5H  =  - . 55 
FROL  =  5000  .  0 
C  =  2  . 5 
C 

C  BASE  ELATE  THICKNESS  LOOP 
C 

111  =  0.375 
DO  100  J=1 , 9 

HR IT E( 10,703)  'BASE  PLATE  THICKNESS,  H1=',H1 
703  F  ORMAT ( A27 , F9 . 3 ) 

HR  I  I E  < 10, x) 

C 

C  VISCOELASTIC  THICKNESS  LOOP 
C 

H2=0 .015 
1)0  150  V  =  1 , 9 
C 

C  (  0113 1  RAINING  LAYER  THICKNESS  LOOP 
C 

H3-0 . 0625 
DO  20  0  L  =  1  ,  9 

IIRI1E(10,702)  'VISCOELASTIC  LAYER  THICKNESS,  H2=’,H2 
UR  I  1  EC  10,*) 

HR  I TE(  1  0,7  02  )  'CONSTRAINING  LAYER  THICKNESS,  H3=’,H3 
7  02  FORMA  1 ( A35 , F6 . 9 ) 

HR  I  1 C( 1 0 , 

HR  I  1 E(  1  0,7  00)  'TEMP' , 'MODE' , ' FCP ’ , 'ETA2', 'ETAS',  'G2' 

7  00  F  ORMA  T ( A3 , 3X , A5 , 9A1 5 ) 

C 

C  PLATE  MODE  LOOP,  READ  MODAL  FREQUENCY  AND  COMPUTE  HAVE  NUMBER 

C 

HTOT  =  Hi  III2MI3 
DO  300  1=1,8 
RE  AIK  1  1  ,  X  )  lirci) 

C 

C  CALCULATE  MODAL  FREQUENCY  AND  HAVE  NUMBER  OF  BASE  PLATE 
C 

SUB  7  =  5QRT((E1*CHT0T*X3)*GC)-'(12.*(1.-NU1XX2)*RH01XH1)) 
KQR-HPC I ) / 5  U  B  7 
T  P  =  HP ( I )/(2 . xpi ) 

C 

C  TEMPERATURE  LOOP 
C 

T  =  30 . 0 

DO  900  K=1 , 15 
C 

C  CALCULATE  PROPERTIES  OT  VISCOELASTIC  FOR  GIVEN  TEMPERATURE  AND 
C  MODE 
C 

501  FR10  =  l  00 10(FP)-U2.*(T-T0))/(  52  5.41-10) 

FR  =  ] 0 . **<  TRIO ) 


onn-r-non'-inooi-nnnn  non  non  oooooooooo  ono 


A-(IR10-lOr.]0(rROL))/C 

SUB] -C*( (SL*  SH>xA  K  SL-SH)*( 1 . -5QRK 1 .+AXX2 . ) ) )/2. 
ETA210  =  LOG10(ET  FROL  HSUBl 


VISCOELASTIC  LOSS  FACTOR 

EIA2=10.xx(ETA210) 
SUB2=2  .  XLOGIOLMROM/HL  ) 
SUB 3  =  1  .  f  (  FROM/  FR  )  x  X|J 
Ml  0  =  L  0G1 0 ( ML )  +  3UB2/SUB3 

VISCOELASTIC  SHEAR  MODULUS 


G2r 1 0 . xx ( Ml  0 ) 


VISCOELASTIC  YOUIIG'S  MODULUS 


E2-G2X2.XC1 . h  MU2 ) 

CALCUI  AT  I OIIS  FOR  DAMPED  PLATE  AT  TEMPERATURE  T  USING  ROS 
KERITI H-UNGAR  EQUATIOtlS 


0-02/ ( E3XH3XH2XK0R) 

U21  =  (HHH2)/2. 

H31=H2*  ( 111  +H3  )/2  . 

Cl  ^Elxlll  x(  1  .  4G)  »G*E3xH5 
D-GXETA2X(E1XH1+E3XH3) 

A I  PIIRE^GXEI  x||]  XE3XMJXUI31  xx2  .  )  x  ( Cl  +  Dx  ET  A2  ) 

AL TUI  11- Gx E 1  XHl  XE3XH3X 11131 XX2. 1XC1XETA2 
BRE  =  E1X||1XE2X||2XH31X(C1  +  DXETA2) 
B1M=E1vII1xE2xU2xH31X(C1xETA2'D) 

5UB8 -  2 . XGXE2XH2XE3XH3XH21XH31XC 
DELRF.-SUB8x((l  .-(ETA2XX2.  ))  +  (DX(2.xETA2))) 
DFLIM-SUBSx(Clx(2.XEIA2)-(Dx(l.-(tTA2xx2.)))) 
:»UBA  =  (  12  ./(C1XX2  .  I  DX»2.  )  )  ^  ( ALPHRE- BRE- DEL  RE ) 
EHCUBEr(El*(Hlxx3.  >  )  *E3*(H3xx3 .  >  +  SUB<t 

MODAL  FREQUENC  <  OF  DAMPED  TLATE 

DETTS  r  RH01  x|II  4  RH02*H24  It  J»RH03 
SUB5=(  EHCUBE>GC)/( 12. x(l . -HU1XX2. )XDEM5) 
MCP-rQRxSQR i ( SUB5 ) 
rcP  =  UCIV(2.>Pl  ) 

COMPARISON  OF  TP  AND  TCP 

IT  ( AB3( I . -rP/TCP) . LE  .  0.1)  THEN 
G010  500 
ELSE 

FP-rcP 
GOTO  501 
ENDIF 

COMPUIE  SYSTEM  LOSS  FACTOR 

00  5UB6  =  L12./(Cl**2.-tn*v2.  ) )  x  ( AL  PHIM- BIM-DEL IM ) 

E  I  AS  -  (  1  ./  EUCLID  E  )  XSU  B6 

PRINT  RESULTS 

HR  IT  E( 1 0, 7  01  )  T  ,  I  ,  FCT , E T A2 , E T AS , G2 
01  FORMAT C5X.F7  .  3 , 2X ,  I  2 , 3X  ,  0  E 1  5  .  <\  ) 

NEXT  TEMPERATURE 

T  -  T  4  5  .  n 
00  CONTINUE 

NEXT  MODE 


HR  I T  E ( 1 0 , x  ) 
300  CONTINUE 


o  n  >—  nnnro  o  o  o 


REHI ND( UNI T  =  1 1 ) 


NEXT  CONSTRAINING  LAYER  THICKNESS 

H  S  =  H  3 ' 0 . 0625 
no  CONTINUE 

NEXT  VISCOELASTIC  LAYER  THICKNESS 

HZ-N2 '0.015 
50  CONTINUE 

NEXT  BASE  fL AT  E  THICKNESS 


Hl=Hl+n . 125 
100  CONTINUE 

Cl  OSE(tlNIT'IO) 
Cl  05E( UNI T  =  1 1  ) 
END 
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APPENDIX  B 


FORTRAN  PROGRAM  USED  TO  COMPUTE  MODAL  LOSS  FACTORS 
FOR  THE  DOUBLE  DAMPING  LAYER  CONFIGURATION  DESIGN 

This  program  uses  the  Ross-Kerwin-Ungar  equations  of  Chapter  3 
to  compute  estimated  modal  frequencies  and  loss  factors  for  the  double 
layer  configuration.  Modal  frequencies  for  an  undamped  plate  are  read 
from  a  data  file  and  estimated  modal  frequencies  and  loss  factors  for 
various  layer  thickness  combinations  in  the  double  layer  configuration  are 
output  to  another  data  file.  Material  properties  of  ISD-112  are  computed 
using  University  of  Dayton  data  and  curve-fitting  equations  to  the  reduced 
frequency  nomogram  [Ref.  9,12].  The  units  used  in  this  program  are 
pounds,  inches,  and  seconds. 
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no  o  n  nnnnnnnnnnnn  a  nr~.  rjoooooooooooonoooonoooonoonnoooonoooonoononn 


PROGRAM  1  MOL YR 


XXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

THIS  PROGRAM  COMPUTES  THE  SYSTEM  I  OSS  PACT  OR  AIID  MODAL  FREQUENCY  OF 
A  DOUBLE  CONSTRAINED  LAYER  V I SCOEL AS T I C AL L Y  DAMPED  PLATE.  THE  LOSS 
FACTORS  ARE  COMPUTED  FOR  A  SPECIFIC  r LA J E/ DAMPING  LAYER 
CONFIGURATION  AND  OVER  A  TEMPERATURE  RANGE  OT  30-100  DEGREES 
FARENIIEI  T  . 

THIS  PROGRAM  APPLIES  TO  A  FRF F - FREE-FREE- FREE  PLATE  AND  THE  UNDAMPED 
MODAL  FREQUENCIES  ARE  DETERMINED  FROM  FINITE  ELEMENT  ANALYSIS  AND 
ARE  READ  INTO  THIS  PROGRAM  FROM  FILE  '  PL  I FRQ ' . 

THE  VISCOELASTIC  MATERIAL  IS  3M  ISD-112.  VISCOELASTIC  MATERIAL 
DATA  IS  FROM  UNIVERSITY  OF  DAYTON  RESEARCH  INSTITUTE 

THE  UNITS  USED  III  THIS  PROGRAM  ARE  ID,  INCH,  SEC,  AND  DEGREES  F. 

LOSS  FACTORS,  DAMPED  PLATE  MODAL  FREQUENCIES,  AND  ISD-112  FROTERT I ES 
ARE  COMPUTED  IN  SUBROUTINE  ' RKU '  FOR  EACH  ll-TH  CONSTRAINED  LAYER 
SYSTEM.  UPON  COMPLETION  OF  EACH  TEMPERATURE  COMPUTATION  THE  SYSTEM 
LOSS  PACT  OR  AND  CORRESPONDING  DAMPED  TLATE  FREQUENCY  ARE  WRITTEN  TO 
TILE  'TITOLYR  DATA' . 

THE  F  01  L  ON  I  III »  COErFICIFIITS  ARC  DEFINED: 

El  =  YOUNG'S  MODULUS  OF  BASF  PLATE  (PSI) 

E  j  =  EQUIVATENT  YOUNG'S  HODIII  US  OF  N-TH  CONSTRAINED  LAYER  SYSTEM 
E3PP.T1  =  YOUNG'S  MODULUS  OT  INF  CONSTRAINING  LAYER  IN  N-TH  LAYER 
ETICUBE  -  EQUIVALENT  STimiEGS  OF  CONSTRAINED  LAYER  SYSTEM 
ETAS  =  SYSTEM  LOSS  FACTOR 

E  t  A  j  -  LOSS  FACTOR  OT  IT-111  COIISTRAINFD  LAYER  SYSTFM 

ElAIFM  --  LOSS  FACTOR  OF  CONSTRAINING  LAYER  IN  THE  N-TH  LAYER  (=0 

TCP  =  FREQUENCY  OF  HIE  DAMPED  PLATE  (HZ) 

TP  =  FREQUENCY  Of  UNDAMPED  PLATE  (HZ) 

GC  =  GRAVITATIONAL  CONSTANT 

til  -  THICKNESS  OF  BASE  TLATE 

H2  =  THICKNESS  OF  1ST  VISCOELASTIC  LAYER 

IT  5  =  THICKNESS  OF  N-IH  CONSTRAINED  LAYER  SYSTEM 

IIA  =  TOIAT  DAMPED  PLATE  THICKNESS 

II 1 PRM  =  THICKNESS  OF  BASE  LAYER  III  N-TH  LAYER 

H2PRM  -  THICKNESS  OF  VEM  IN  N-TH  LAYER 

H3PRM  -  THICKNESS  OF  CONS  1  R  A I N I  TIC-  LAYER  IN  N-TH  LAYER 

KQR  =  HAVE  NUMBER  OF  UNDAMPED  PI  ATE 

IT  n 1  =  POISSON'S  RATIO  OF  BASE  FI  A I F  AND  CONSTRAINING  LAYERS 

IIM2  -  rOISSON'S  RATIO  OT  V I SCOEL AS  I  I C  MATERIAL 

HN2PRM  =  POISSON'S  RATIO  OF  VEM  IN  N-TH  LAYER 

R HOI  =  DENSITY  OT  BASE  fM  ATE  AND  CONSTRAINING  LAYERS 

RIIU2  -  DENSITY  OT  VISCOELAS1IC  MATERIAL 

RHUS  PIT  =  DCNSITY  OF  VEM  III  NTH  LAYFR 

RII0  3  =  DENSITY  OF  N-IH  LAYER  (=  RII01) 

RH03PM  =  DENSITY  OF  ll-TH  LAYER  CONS  TRAINING  LAYER 
1  -  T  EMPERA I URF  VARIABIE 

HP  =  FREQUENCY  OF  UNDAMPED  TLA1E  (RAD/SEC) 

C,F.IFROI  .  r  ROL  ,  F  ROM .  ML  ,  ITROM ,  II ,  SH ,  SI  ,T0  =  COErriCIEIHS  FOR  THE 
REDUCED  FREQUENCY  NOMOGRAM  EQUATIONS 

XXX*X***XX**XXXXXXXXX*XXX¥XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

REAI  C,  FI  .  E  3  ,  E  3  T  F.  M ,  F.MCUBE ,  FT  AS  ,  FT  A  3 ,  r  T  A3rM ,  ET  FROL  ,  TCP,  FP,  FROL 
R  E  A  1  FROM. GC, 111  ,  M2  ,  H  3  .  HA  ,  H 1  TRII ,  II2PRM  ,  H3PRM  .  KQR  ,  ML  , ITROM, H,NU1 
REAI  1 1 1 1 2 , 1 1 U  2  P  r  1 ,  P  I  ,  R 1 1 U I  ,  RII02,RII02rn,  SH,  SL  ,  SUB1  ,  T  ,  T  0  ,  IIP 
REAL  P.H03 ,  RH03PM 

DIMENSION  l  IT  (  17  ) 

PI  -'t .  0*A  I  AIK  1  .  0  ) 

OPLN(HNIT-]0,FIIE-'ri TFRQ' , STATUS-' 01 D' ) 

OP  EN(  UN  I  T  -  1  1  ,  T  1  L  t.  -  '  IIIULYR  '  , STATUS-' OLD'  ) 

ASSIGN  LAYER  THICKNESSES 


8  5 


non  non  non  n  n  n  -j  o  o  o  n  n  n  o  non 


111  -0  .  5 

II?  -  0  .  0 4  5 

llll'Rlt=0  .09375 

M2rRH=  0 .015 

H3PRM=0  .  09575 

115  =  111  I'RM  l  II 2 PR IH  II5PRM 

II ■  i  ^  It  1  i  H2  <  H3 

ASSIGN  MATERIAL  C0NSIAN1S 

E  1  -  1  .  0  E  7 
RIIU1  =0 . 096  8 
111)1  =  0 . 33 
RH02  =  0 . 035 
11112  =  0  .  59 
RHU2PM  =  0 .0  35 
H 1 1 2  R  IT  -  0  .  5 
RH03-0 .0968 
E5PRH-1 . 0E7 
RIIUSPM  =  0  .  0968 
ETA  3PM  -  0 . 0 
El  A3  =  0.0 

GC-  386  .  0 

DEEIIIE  VISCOELASTIC  COHSTAIITS  TOR  NOMOGRAM  EQUATIONS 

10=  106.0 
T  ROM- 2 . 0E6 
MR(.'M  =  6 38 . 96 
11  =  0  .27  5 
ML =8  .  7 
L  I  f  R  O I.  -1.08 
SI -0.65 
Sll  -  -0  .  55 
FRUI =5000.0 
C  =  2 . 5 


IJRITE  PLAIE  CHARACTERISTICS 


00 
0  1 


iiri  irm  ,700) 

HR  II  EC  1  1  ,  7  00  ) 
IIRI  If  (  1  1  .  7  00  ) 
IIRI  I  U  1  1  .  7  00  ) 
I IR I  II  C  1  1  ,  7  0  0  ) 
IIRI  IE(  1  1  ,* ) 
mu i m  i .  * ) 


111  r  ’  ,  H 1 
HE  -  •  .112 
HI  PRIME  = 
II2PRIME  = 
II5PRIME  = 


ill  r  rm 

II2PRM 
H  51' RM 


I  ORMA I ( A  1 2 , E6 . 5  ) 

IIRI  I  !  C  1  I  ,  701  )  '  I  EMP'  ,  'MODE'  ,  'TCP 
EOKMA I ( AS , 3X, A 5 , 2A1 5 ) 

hr i reel l , >  ) 


'ETAS' 


PLAIE  MODE  LOOT 


DO  100  1-1,17 
READ!  10,  x)  IIP (  I  ) 


CALCUL A I E  MODAL  FREQUENCY  AND  HAVE  NUMBER  Or  UNDAMPED  PLATE 

sum  -  SORT  (  (  E1*(IK,XX3)XGC)/U2  ,X(  1  .  -NU1  x*2  )  XRH01  *116  )  ) 
KQR-IIP  (  I  ) / S U B 1 

r  p  =  iip(  i )/( 2 .  xri ) 


T  EMP  ERA  I URE  LOOP 
I  -  50  .  0 

DO  20  0  K-J  .  15 

COMPUIE  SYS  I  EM  LOSS  FACTOR  AND  S1IPFNESS  FOR  N-TH  LAYER 

CAL  L  RKU(  HI  PRM,  II2PRM,  H3PRM,  El  ,  RHO 1  ,  HU  1  ,  NU2PM ,  RH02TM ,  E3TRM,  RH03PII, 
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cr i  asi  m,  t , rr, kqr, et as, eiicube, rcr,  to,  prom, mrom, n,ml , etfrol , sl ,sh, 

C  P  R  U  l  ,  C  ) 

c 

C  CQNYER 1  RCSUL  T  3  PROtl  H-Ilt  LAYER  CALCULAriOIl  10  T01AL  PLAT  E 

C 

E  T  A  3  -  0.0 

pp  =  lire  i )/  c  2 .  *n  ) 

c  pp-rr.r 

E  3  =  EIICUBE/  (  H3*x3  ) 

C 

C  COMPUTE  SYSTEM  LOSS  PACTOR  AMD  FREQUENCY  FOR  TOTAL  PLATE 
C 

CALL  RKU(  111  ,112, 113,  El  ,  RIIOl  ,  IIN1  ,  NU2  ,  RH02  ,  E  3 ,  R II 0  3  ,  ETA3.T,  FP,KQR,  ETAS, 
CEIICUBE  ,  TCP  ,  T  0  ,  r  ROM  ,  MROM ,  N  ,  ML  ,EIFR0L,SL,5II,PR0L,C) 

C 

C  nil  NT  RESULTS 

URITE(11,702)  T,l, Per, ETAS 
7  02  FORMA  T  (  SX  ,  F 7  .  3 , 2X ,  I  2 , 3X ,  2E1 5  .  <t ) 

C 

C  NEXT  TEMPERATURE 
C 

T  -  1  i  5  .  0 

roo  con  r i hue 

c 

C  NEXT  MOPE 

c 

HR  I T  E(  1  1 ,  X ) 

ion  con  t  j  hue 

c 

CLOSECUM! 1-11) 

Cl  USE (UNI  1=10) 

END 

c 

o  ***'»*'(XX»vxXXXX*XXXXy*xXXX»XXXXXXX»'*XXXXXXX*XXXXXXXXXXXXXXXXXXXXXXXXXX 
C  yvyvimxXvXX'KXXXXXXXXXXXXXXXXXXXXXXXXX XXXXKXKXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

SUBROUTINE  RKU  ( H 1 ,  H?  ,  H  3 ,  E 1 ,  RHO 1  .  Mill ,  NU2 ,  RH02 ,  E3  ,  RU03 ,  E 1 A  3 ,  T  ,  FT  , 

C  KQR, ETAS, EIICUBE, FCP, TO, FROM,  MROM, 11,  ML, ETFROL, SL,SH, FROL, C) 

C 

C  XXXVXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

C  THIS  SUP RON  1 1 NE  CALCULATES  VISCUEIASTIC  PROPERTIES  BASED  ON  THE 
C  UNIVF.F5I1F  OF  DAYTON  REDUCED  FREQUENCY  NOMOGRAM  EQUATIONS,  AND  THEN 

C  CALCULATES  PLATT  S I  1 EFNESSES  AND  LOSS  FACIORS  BASED  ON  THE  ROSS- 

C  KtRITI  N-UNGAR  EQUATIONS. 

C 

C  IMF  FULL  (THING  ADDITIONAL  VARIABLES  ARE  DEFINED  TOR  USE  IN  THIS 

C  SNBRUU  HUE: 

c 

C  A  =  COFF  T ICI ENT  TOR  NOMOGRAM  EQUATIONS 

C  ALPII1M, ALTHRE  =  IMAGINARY  AND  REAL  COMPONENTS  OF  COEFFICIENT 

C  ALPHA  IN  THE  RKU  EQUATIONS 

C  BIP1, ORE  =  IMAGINARY  AND  REAL  COMPONENTS  OF  COEFFICIENT  'B'  IN 

C  I  HE  RKU  E QUA  I  I  0115 

c  ci.d  -  coet ricir.tits  fop.  rku  equations 

C  DELHI. DELRE  =  IMAGINARY  AND  REAL  COMPONENTS  OF  COEFFICIENT 

C  DELIA  IN  THE  RKU  FQUAI10NS 

C  DENS  =  C  0  f  1 B I  N  A  1  ION  OF  MAIF.RIAI  DENSITIES  U5ED  TO  COMPUTE  THE 

C  FREQUENCY  OF  THE  DAMPED  PLATE  AND  DAMPING  LAYERS 

C  El ICUB E  =  FQUIVALRII  STIFFNESS  OF  DAMPED  PLATE  AS  COMPUTED 

C  USING  RKU  EQUATIONS 

C  F  T  A  2  =  LOSS  r  A  C I  0  R  or  VEM  COMPUTED  IN  NOMOGRAM  EQUATIONS 

C  El  AD  10  =  L  0  G 1  0  (  F.  I  A  2  ) 

C  E I  AS  -  SYSTEM  LOSS  FACTOR  COMPUTED  BY  RKU  EQUATIONS 

C  TCP  -  MODAL  IRE OUFNCY  OE  DAMPED  PLATE  (IIP) 

C  IP  =  MODAL  FREQUENCY  OE  UNDAMPED  PLAIE  (HZ) 

C  I R  =  REDUCED  FREQUENCY  OF  VEM 

C  FRIO  -  L  OG1  l)(  FR  ' 

C  G2  -  SHEAR  MODULUS  OF  VEM 

C  1121. NT]  -  RKU  EQUAIION  COET  F  I  Cl  Fll  T  S 

C  M 1 0  =  LOGIOCG2)  AS  COMPUIED  BX  NOMOGRAM  EQUATIONS 
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non  non  n  non  non  nnn  noo  mnnn  o  on  on 


lice  -  FREQUENCY  or  DAMPED  PLATE  CRAD/SEC) 

XXXKXVXXVXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXX 

REAL  A.  A I  rniM, A L PURE, DIM, BRE.C.Cl , D, DELHI, DEI  RE , DENS, El , E2, E3 
RLAl  EMCUBE,  FT  A?,  Et  A3,  ET  121  0,  E  I  AS,  E  I  FROL  ,  FCP,  TT,  TP.,  FRIO,  FROl ,  FROM 
REAL  G . G2 , GC , H 1  ,112,(13,1121  ,  H3  1  ,  KQR ,  111  0  ,  ML  ,  HKOH ,  N ,  NU1  ,IIU2,  PI  ,RI10I 
RE  A  |  RII02,  RIIOJ,  311,  SI  ,  SUB  1  ,  SU B2 ,  SUB  3 ,  SUBA ,  SUB5 , 5UB6 , SUD7  ,  SUB8 
RLAL  SUD9. SUB10, SUB1 1 , T , T 0 ,HCP 

PI  -A  .  0*ATAtl(  1.0) 

GC=  386  .  0 

CALCULATE  PROPERTIES  OF  VEH  TOR  GIVEN  TEMPERATURE  AMD  MODE 

1  rR10-LOG10(FP)-(12.x(T-T0))/(525.-H-TO) 

FRHO.xx(iRlO) 

A  =  (i'R10-tUG10(FROL)  )/C 

SUD1 =C*( ( SL+SH)*A*( SL-SH)x( 1 ,-SQRTL 1 . +AXX2) ) )/2 . 

ETA210-LOG10C  ETFP.OL  I  +  SUD1 

VISCOELASTIC  LOSS  FACTOR 

E  T  A  2  =  1  0  . x  x  ( ETA210) 

VISCOELASTIC  SHEAR  MODULUS 

SUB2  =  2 . X  L  0G1 0 ( MROM/ML ) 

SUB3  =  1  .  +  C  FROM/  FR  )  X x(T 
Ml 0=LOG10(ML )+SUD2/SUB3 
G2  =  1 0  . xx ( Ml  0 ) 

VISCOELASTIC  YOUNG'S  MODULUS 
E2  rG2x 2  .  x (  1  .  IIU2  ) 

CALCULATIONS  FOR  DAMPED  PLATE  USING  RKU  EQUATIONS 

G-G2/I  E3XH3XH2*ILQR) 

H21=CHl+H2)/2. 

H31=H2+(Hl+H3)/2. 

Cl=ElxHlx(l . +G)+GXE3XH3X( 1 .-ETA2XETA3) 

D=GXE1XH1 XETA2«G*E3*H3X(ETA24ETA3) 

SUBA=G*E1*H1xe3xH3X(H31xx2) 

ALPMRE=SUBAX(C1X(1.-ETA2*ETA3)+DX(ETA2+ETA3)) 

ALPHIM  =  S’JBAX(C1X(  ETA2IETA3)) 

SUB5-E1XH1XE2XH2XH31 

BRE-SUB5x(Cl+DxETA2) 

BIM  =  SUB5X(C1XE.TA2-D) 

SUD6=2.XG*E2*H2*E3*H3*H21xH31 
SUB  7  =  1 .-2.*ETA2*ETA3-(ETA2xx2) 

SUBS =2 . XETA2+ETA3-C  E T A2xx2 ) x E T A3 
DELRE-SUB6x(ClxSUB7  +  DxSUB8  ) 

DEL IM-SUB6XCC1XSUB8-DXSUB7 ) 

SUB9=(12./<Clx*2*D«*2) ) x ( AL PHRE- BRE- DE L RE ) 
ENCUDE=Elx(Mlxx3)tE3x(H3xx3)+SUB9 

MODAL  FREQUENCY  OF  DAMPED  TLATE 

DENS- RHOl XH14RH02*H2+RH03xH3 

SUBI0  =  (  EHCUBExGC  )^(  1  2  .  x ( IIU1  XX2  ) x DENS ) 

WCP  =  KQRxSQRi(SUBlO) 

FCP  =  HCP/(  2  .  xn  ) 

COMPARISON  OF  FP  AND  TCP 

I rc absc i . -FP/rcr )  .le.  o.io)  then 

GOTO  500 

El  SE 

Fp=rcp 

GOTO  501 
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nonn 


lifIDI  F 


COMTE)  7 E  SYSTEM  LOSS  rACTOR 

SUBll  =  (12./'(Clx*2<D**2))*(ALPHIM-BIM-DELIM) 
ETAS=( 1 . /EHCURE)*( (E3>H3*ETA3)  +  SUBn) 

C 

END 
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APPENDIX  C 


DESIGN  DRAWINGS  FOR  THE  MACHINING 
OF  THE  FLOATING  ELEMENT  AND 
POCKET  PLATE  CONFIGURATIONS 


The  drawings  shown  in  Figures  C.l  and  C.2  were  used  to  machine 
the  pocket  plate  and  floating  element  plate  used  in  the  experiments.  These 
drawings  are  included  to  show  the  relation  of  the  pocket  for  the  ISD— 112  to 
the  cover  plate  and  how  the  viscoelastic  was  protected  from  the  heat  of 
welding. 
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Figure  C.2.  Design  drawing  of  the  floating  element  configuration. 
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APPENDIX  D 


REPRESENTATIVE  MSC/NASTRAN  DATA  DECK  FOR 
THE  DAMPING  CONFIGURATIONS 

This  data  deck  was  used  to  compute  the  modal  frequency  response  of 
the  single  layer  damping  configuration  and  is  a  representative  sample  of 
the  NASTRAN  decks  used  for  the  other  finite  element  models.  The  values 
in  the  damping  table  are  from  a  curve  fit  to  the  modal  loss  factors  estimated 
from  the  modal  strain  energy  method.  Since  the  data  deck  for  the  normal 
mode  and  modal  strain  energy  extraction  is  virtually  identical  to  this  deck, 
the  Case  Control  deck  commands  for  the  normal  mode  analysis  are 
included,  but  are  commented  out. 

The  OUTPUT  request  provides  data  for  an  x-y  plot  of  the  modal 
frequency  response. 

The  units  used  in  this  deck  are  pounds,  inches,  and  seconds. 
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Id  single, ■ f r 
sol  30 
$ 

5  iimiiiiimiiimmiiiiiiiiiiiiiiiiiiiiiiiiHiiiiiiiiiiimimmil 

$  THIS  DECK  IS  TO  COMPUTE  THE  HODAI,  FREQUENCY  RESPONSE  OF  A  15i45 

$  INCH  ALUMINUM  PLATE  WITH  A  CONSTRAINED  VISCOELASTIC  DAHP1NG  LAYER. 

$  THE  HODEL  HAS  252  ELEKEHTS  WITH  84  ELEHENTS  III  EACH  LAYER.  THE  QUAD 
$  ELEMENTS  ARE  OFFSET  FROM  THE  HEY  ELEMEMTS  AS  SUGGESTED  BY  THE 
$  LITERATURE. 

$ 

$  THE  MATERIAL  PROPERTIES  OF  ISD-112  ARE  FROH  THE  3H  CORPORATION 

$ 

$  THE  MODAL  LOSS  FACTORS  II  THE  DAMPING  TABLO  ARE  FROM  A  CURVE-FIT 

$  TO  THE  SET  OF  MODAL  LOSS  FACTORS  COMPUTED  FROH  THE  MODAL  STRAII  ENERGY 
$  METHOD. 

$ 

$  THE  UNITS  USED  II  THIS  DECK  ARE  POUIDS, INCHES,  AID  SECONDS 

J  1  tilt llllll till  tilt lllllllElllltlllltlll <<1111 < Hill 1 1IIIII<I1<I111<1I< 

tl«e  3000 
cend 

title  =  MODAL  FREQUENCY  RESPOISE  /  200  HZ  /  3H 

■ethod  =  1 

spe  =  1 

dload  -  10 

frequency  =  10 

sdaiping  =  101 

set  111  -  95  , 

svector  =  all 

acceleratlon(plot, phase)  =  111 
output { xyplot ) 
xyprlnt  acce  /  95(t3) 

$ 

$  iiiiiiiitiiiiii 

$  THE  FOLLOWING  LINES  ARE  THE  CASE  CONTROL  DECK  CARDS  FOR  THE 
$  NORMAL  MODE  EXTRACTION  AND  STRAII  EIERGY  REQUEST 
$  iethod  =  1 

$  3pC  =  1 
$  set  10  =  all 
$  set  11  =  169, thru, 252 
$  ese  =  11 
$ 

$  I<<<I<<IIII1I<< 

BEGIN  BULK 


$  TITLE  =  SINGLE  LAYER 

WITH  QUAD  OFFSET 

$  DATA  DECK 

PRODUCED  BY 

PATNAS  VERSION  2.0:  24-IOV 

GRID 

1 

45.0000  15.0000  0.51500 

GRID 

2 

45.0000  15.0000  0.50000 

GRID 

3 

45.0000  12.5000  0.51500 

GRID 

4 

45.0000  12.5000  0.50000 

GRID 

5 

45.0000  10.0000  0.51500 

GRID 

6 

45.0000  10.0000  0.50000 
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GRID 

7 

45.0000 

7.50000 

0.51500 

GRID 

8 

45.0000 

7.50000 

0.50000 

GRID 

9 

45.0000 

5.00000 

0.51500 

GRID 

10 

45.0000 

5.00000 

0.50000 

GRID 

11 

45.0000 

2.50000 

0.51500 

GRID 

12 

45.0000 

2.50000 

0.50000 

GRID 

13 

45.0000 

0. 

0.51500 

GRID 

14 

45.0000 

0. 

0.50000 

GRID 

IS 

41.7857 

15.0000 

0.51500 

GRID 

16 

41.7857 

15.0000 

0.50000 

GRID 

17 

41.7857 

12.5000 

0.51500 

GRID 

18 

41.7857 

12.5000 

0.59000 

GRID 

19 

41.7857 

10.0000 

0.51500 

GRID 

20 

41.7857 

10.0000 

0.50000 

GRID 

21 

41.7857 

7.50000 

0.51500 

GRID 

22 

41.7857 

7.50000 

0.50000 

GRID 

23 

41.7857 

5.00000 

0.51500 

GRID 

24 

41.7B57 

5.00000 

0.50000 

GRID 

25 

41.7857 

2.50000 

0.51500 

GRID 

26 

41.7857 

2.50000 

0.50000 

GRID 

27 

41.7857 

0. 

0.51500 

GRID 

28 

41.7857 

0. 

0.50000 

GRID 

29 

38.5714 

15.0000 

0.51500 

GRID 

30 

38.5714 

15.0000 

0.50000 

GRID 

31 

38.5714 

12.5000 

0.51500 

GRID 

32 

38.5714 

12.5000 

0.50000 

GRID 

33  . 

38.5714 

10.0000 

0.51500 

GRID 

34 

38.5714 

10.0000 

0.50000 

GRID 

35 

38.5714 

7.50000 

0.51500 

GRID 

36 

36.5714 

7.50000 

0.50000 

GRID 

37 

38.5714 

5.00000 

0.51500 

GRID 

38 

38.5714 

5.00000 

0.50000 

GRID 

39 

38.5714 

2.50000 

0.51500 

GRID 

40 

38.5714 

2.50000 

0.50000 

GRID 

41 

38.5714 

0. 

0.51500 

GRID 

42 

38.5714 

0. 

0.50000 

GRID 

43 

35.3571 

15.0000 

0.51500 

GRID 

44 

35.3571 

15.0000 

0.50000 

GRID 

45 

35.3571 

12.5000 

0.51500 

GRID 

46 

35.3571 

12.5000 

0.50000 

GRID 

47 

35.3571 

10.0000 

0.51500 

GRID 

48 

35.3571 

10.0000 

0.50000 

GRID 

49 

35.3571 

7.50000 

0.51500 

GRID 

50 

35.3571 

7.50000 

0.50000 

GRID 

51 

35.3571 

5.00000 

0.51500 

GRID 

52 

35.3571 

5.00000 

0.50000 

GRID 

53 

35.3571 

2.50000 

0.51500 

GRID 

54 

35.3571 

2.50000 

0.50000 

GRID 

55 

35.3571 

0. 

0.51500 

GRID 

56 

35.3571 

0. 

C.SCOOO 

GRID 

57 

32. i 429 

15.0000 

0.51500 
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GRID 

be 

32.1429 

15.0000 

0.50000 

GRID 

59 

32.1429 

12.5000 

0.51500 

GRID 

60 

32.1429 

12.5000 

0.50000 

GRID 

61 

32.1429 

10.0000 

0.5)500 

GRID 

62 

32.1429 

10.0000 

0.50000 

GRID 

63 

32.1429 

7.50000 

0.51500 

GRID 

61 

32.1429 

7.50000 

0.50000 

GRID 

65 

32.1429 

5.00000 

0.51500 

GRID 

66 

32.1429 

5.00000 

0.50000 

GRID 

67 

32.1429 

2.50000 

0.51500 

GRID 

68 

32.1429 

2.50000 

0.50000 

GRID 

69 

32.1429 

0. 

0.51500 

GRID 

70 

32.1429 

0. 

0.50000 

GRID 

71 

28.9286 

15.0000 

0.51500 

GRID 

72 

28.9286 

15.0000 

0.50000 

GRID 

73 

28.9286 

12.5000 

0.51500 

GRID 

74 

28.9286 

12.5000 

0.50000 

GRID 

75 

28.9286 

10.0000 

0.5i500 

GRID 

76 

28.9286 

10.0000 

0.50000 

GRID 

77 

28.9286 

7.50000 

0.51500 

GRID 

78 

28.9286 

7.50000 

0.50000 

GRID 

79 

28.9286 

5.00000 

0.51500 

GRID 

80 

28.9286 

5.00000 

0.50000 

GRID 

81 

28.9286 

2.50000 

0.51500 

GRID 

82 

28.9286 

2.50000 

0.50000 

GRID 

03 

28.9286 

0. 

0.51500 

GRID 

84 

28.9286 

0. 

0.50000 

GRID 

85 

25.7143 

15.0000 

0.51500 

GRID 

86 

25.7143 

15.0000 

0.50000 

GRID 

87 

25.7143 

12.5000 

0.51500 

GRID 

88 

25.7143 

12.5000 

0.50000 

GRID 

89 

25.7143 

10.0000 

0.51500 

GRID 

90 

25.7143 

10.0000 

0.50000 

GRID 

91 

25.7143 

7.50000 

0.51500 

GRID 

92 

25.7143 

7.50000 

0.50000 

GRID 

93 

25.7143 

5.00000 

0.51500 

GRID 

94 

25.7143 

5.00000 

0.50000 

GRID 

95 

25.7143 

2.50000 

0.51500 

GRID 

96 

25.7143 

2.50000 

0.50000 

GRID 

97 

25.7143 

0. 

0.51500 

GRID 

98 

25.7143 

0. 

0.50000 

GRID 

99 

22.5000 

15.0000 

0.51500 

G81D 

100 

22.5000 

15.0000 

0.50000 

GRID 

101 

22.5000 

12.5000 

0.51500 

GRID 

102 

22.5000 

12.5000 

0.50000 

GRID 

103 

22.5000 

10.0000 

0.51500 

GRID 

104 

22.5000 

10.0000 

0.50000 

GRID 

105 

22.5000 

7.50000 

0.51500 

GRID 

106 

22.5000 

7.50000 

0.50000 

GRID 

107 

22.5000 

5.00000 

0.51500 

GRID 

108 

22.5000 

5.00000 

0.50000 
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22.5000 
22.5000 
22.5000 
22.5000 
19.2857 
19.2857 
19.2857 
19.2857 
19.2057 
19.2857 
19.2857 
19.2857 
19.2857 
19.2857 
19.2857 
19.2857 
19.2857 
19.2857 
16.0714 
16.0714 
16.0714 
16.0714 
16.0714 
16.0714 
16.0714 
16.0714 
16.0714 
16.0714 
16.0714  ; 
16.0714  ; 
16.0714  l 
16.0714  I 
12.6571 
12.8571  . 
12.8571  ! 
12.8571  I 
12.8571  1 
12.8571  j 
12.8571  1 
12.8571  ( 
12.8571  5 
12.8571  5 
12.8571  2 
12.8571  2 
12.0571  0 
12.8571  0 
9.64286  1 
9.64286  1 
6.42057  1 
6.42857  1 
3.21429  1 


•0  2.50000  0 
0  2.50000  0 
0  0.  0 
0  0.  0 
7  15.0000  0 
7  15.0000  0 
7  12.5000  0 
7  12.5000  0, 
7  10.0000  0 
7  10.0000  0, 
7  7.50000  0. 
7  7.50000  0, 
7  5.00000  0. 
7  5.00000  0, 
7  2.50000  0. 
7  2.50000  0. 
7  0.  0. 

7  0.  0. 

I  15.0000  0. 
I  15.0000  0. 
I  12.5000  0. 
1  12.5000  0. 
I  10.0000  0. 
I  10.0000  0. 
I  7.50000  0. 

I  7.50000  0. 

I  5.00000  0. 
5.00000  0. 

2.50000  0. 

2.50000  0. 
0.  0.! 
0.  0.! 
15.0000  0.! 
15.0000  0.! 
12.5000  0.! 
12.5000  0.! 
10.0000  0.! 
10.0000  o.; 

7.50000  0.5 

7.50000  0.5 
5.00000  0.5 
5.00000  0.5 

2.50000  0.5 

2.50000  0.5 
0.  0.5 

0.  0.5 

15.0000  0.5 
15.0000  0.5 
15.0000  0.5 
15.0000  0.5 
15.0000  0.5 


0.51500 
0.50000 
0.51500 
0. 50000 
0.51500 
0.50000 
0.51500 
0.50000 
0.51500 
0.50000 
0.51500 
0.50000 
0.51500 
0.50000 
0.51500 
0.50000 
0.51500 
0.50000 
0.51500 
0.50000 
0.51500 
0.50000 
0.51500 
0.50000 
0.51500 
1.50000 
0.51500 
0.50000 
• . 51500 
.50000 
.51500 
.50000 
.51500 
.50000 
.51500 
.50000 
.51500 
.50000 
.51500 
.50000 
.51500 
.50000 
.51500 
.50000 
51500 
.50000 
51500 
50000 
51500 
50000 
51500 


GRID 

160 

3. 21429  15.0000  0.50000 

GRID 

161 

0.  15.0000  0.51500 

GRID 

162 

0.  15.0000  0.50000 

GRID 

163 

9.64286  12.5000  0.51500 

GRID 

164 

9.64286  12.5000  0.50000 

GRID 

165 

9.64206  10.0000  0.51500 

GRID 

166 

9.64286  10.0000  0.50000 

GRID 

167 

9.64206  7.50000  0.51500 

GRID 

166 

9.64286  7.50000  0.50000 

GRID 

169 

9.64286  5.00000  0.51500 

GRID 

no 

9.64286  5.00000  0.50000 

GRID 

171 

9.64286  2.50000  0.51500 

GRID 

172 

9.64286  2.50000  0.50000 

GRID 

173 

9.64286  0.  0.51500 

GRID 

174 

9.64286  0.  0.50000 

GRID 

175 

6.42857  12.5000  0.51500 

GRID 

176 

6.42857  12.5000  0.50000 

GRID 

177 

3.21429  12.5000  0.51500 

GRID 

178 

3.21429  12.5000  0.50000 

GRID 

179 

0.  12.5000  0.51500 

GRID 

160 

0.  12.5000  0.50000 

GRID 

181 

6.42857  10.0000  0.51500 

GRID 

182 

6.42857  10.0000  0.50000 

GRID 

183 

6.42857  7.50000  0.51500 

GRID 

184 

6.42857  7.50000  0.50000 

GRID 

185 

6.42857  5.00000  0.51500 

GRID 

186 

6.42857  2.50000  0.51500 

GRID 

187 

6.42857  5.00000  0.50000 

GRID 

188 

6.42857  2.50000  0.50000 

GRID 

189 

6.42657  0.  0.51500 

GRID 

190 

6.42857  0.  0.50000 

GRID 

191 

3.21429  10.0000  0.51500 

GRID 

192 

3.21429  10.0000  0.50000 

GRID 

193 

0.  10.0000  0.51500 

GRID 

194 

0.  10.0000  0.50000 

GRID 

195 

3.21429  0.  0.51500 

GRID 

196 

3.2M29  0.  0.50000 

GRID 

J97 

0.  0.  0.51500 

GRID 

198 

0.  0.  0.50000 

GRID 

199 

3.214P9  7.50000  0.51500 

GRID 

200 

3.21429  5.00000  0.51500 

GRID 

201 

3.21429  2.50000  0.51500 

GRID 

202 

3.21429  7.50000  0.50000 

GRID 

203 

3.21429  5.00000  0.50000 

GRID 

204 

3.21429  2.50000  0.50000 

GRID 

205 

0.  7.50000  0.51500 

GRID 

206 

0.  5.00000  0.51500 
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